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Abstract. 

We study the statistics of records of a one-dimensional random walk of n steps, 
starting from the origin, and in presence of a constant bias c. At each time-step 
the walker makes a random jump of length r\ drawn from a continuous distribution 
/(?/) which is symmetric around a constant drift c. We focus in particular on the 
case were f{rf) is a symmetric stable law with a Levy index < /i < 2. The record 
statistics depends crucially on the persistence probability which, as we show here, 
exhibits different behaviors depending on the sign of c and the value of the parameter 
/i. Hence, in the limit of a large number of steps n, the record statistics is sensitive 
to these parameters (c and fi) of the jump distribution. We compute the asymptotic 
mean record number (R n ) after n steps as well as its full distribution P(R,n). Wc 
also compute the statistics of the ages of the longest and the shortest lasting record. 
Our exact computations show the existence of five distinct regions in the (c, < /i < 2) 
strip where these quantities display qualitatively different behaviors. We also present 
numerical simulation results that verify our analytical predictions. 
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1. Introduction 

The statistical properties of record-breaking events in stochastic processes have been 
a popular subject of research in recent years. The theory of records has found 
many interesting applications. Record events are very important in sports (TJ [2] and 
climatology [31 HJ El E] , but have also been found relevant in biology [7] , in the theory 
of spin-glasses [H [9] and in models of growing networks [10]. Also in finance, record- 
breaking events, e.g., when the price of a stock breaks its previous records, can lead 
to increased financial activities [TT] IT2] . In all of these fields researchers have recently 
made progress in understanding and modeling the statistics of records by comparing the 
records in observational data with various kinds of stochastic processes. In this context it 
has become increasingly important to improve our understanding of the record statistics 
of elementary stochastic processes. In this paper we focus on one such elementary 
stochastic process namely a random walk in presence of a constant bias. We show that 
even for such a simple process, its record statistics is considerably nontrivial and rich. 

In general, one is interested in the record events of a discrete-time series of random 
variables (RV's) Xo, Xi, x n . An (upper) record is an entry Xk, which exceeds all 
previous entries: Xk > max (xq, x%, Xk-i). Until the end of the last century record 
statistics was fully understood only in the case when the entries of the time series are 
independent and identically distributed (i.i.d.) RV's (see for instance [131 [TH [15] ) . F° r 
i.i.d. RV's from a continuous distribution p(x) the probability r n of a record in the n-th 
time step is given by p3] 

r n := Prob [x n > max (x , Xt, x n _i)] = ^qjy , (1) 

which is universal, i.e., independent of the parent distribution p(x). This universality 
follows simply from the isotropy in ordering, i.e., any one of the (n+ 1) entries is equally 
probable to be a record. Let R n denote the total number of records up to step n. The 
mean record number is then simply (R n ) = J2m=o Tm ' wmcn grows asymptotically as 
~ In n for large n. 

Due to the numerous applications of the theory of records it became interesting to 
consider more general models. There has been a lot of interest in the record statistics 
of RV's which are uncorrelated but not identical anymore. For instance Ballerini et al. 
considered uncorrelated RV's with a linear drift [IE]. More recently Franke et al. studied 
the same problem as well and found numerous new results [TTl [T8l [T9] by also considering 
the correlations between individual record events. This model was then successfully 
applied to the statistics of temperature records in the context of global warming [5] . In 
2006 Krug studied the statistics of records of uncorrelated RV's with a time- increasing 
standard deviation, a model with important biological implications [20] . 

Another important issue is the study of record statistics for correlated random 
variables. For weak correlations, with a finite correlation time, one would expect that 
the record statistics for a large sequence to be asymptotically similar to the uncorrelated 
case. This is no longer true when there are strong correlations between the entries. 



Record statistics and persistence for a random walk with a drift 



3 



Perhaps, one of the simplest and most natural time series with strong correlations 
between its entries corresponds to the positions of a one dimensional random walk |21j . 
Despite the striking importance and abundance of random walk in various areas of 
research, the record statistics of a single, discrete-time random walk with a symmetric 
jump distribution was not computed and understood until only a few years ago. In 2008, 
Majumdar and Ziff [22] computed exactly the record statistics of a one dimensional 
symmetric random walk model and showed that the record rate of such a process is 
completely universal for any continuous and symmetric jump distribution, thanks to 
the so called Sparre Andersen theorem [23] . They considered a time series of RV's x m 
given by: 

•^m •Em—l Vmi (2) 

where rj^s are i.i.d. RV's drawn from a symmetric and continuous jump distribution 
/ (n) (it includes even Levy flights where f(rj) ~ l/|r^ +1 with < fj, < 2). Then, the 
record rate r n for such a process is given by the universal formula [22] 

r n = ( 2n V 2 " ^> -L , (3) 
V n J Wirn 



independently of the jump distribution f(rj). They also computed exactly the mean 
record number (Rn) and even its full distribution [22]. In addition, there exists nice 
connection between the record statistics and the extreme value statistics for the one 
dimensional symmetric jump processes and many universal results can be subsequently 
derived using the Sparre Andersen theorem (see [23] for a review). 

Following Ref. [22], there has been considerable interests in generalising them to 
more general set of strongly correlated stochastic processes. For instance, Sabhapandit 
discussed symmetric random walks with a random, possibly heavy tailed, waiting time 
between the individual jumps (the so called Continuous Time Random Walk model) [25] . 
Recently the present authors considered the record statistics of an ensemble of N 
independent and symmetric random walks JT2] . There, in contrast to the case of a 
single random walker, the record statistics of N Levy flights with a heavy-tailed jump 
distribution was found to be different from the one of N Gaussian random walkers with 
a jump distribution that has a finite second moment. 

Another important generalization is to consider a single one dimensional random 
walker but with asymmetric jump distribution, for instance, in presence of a constant 
bias c. First steps towards this generalization were taken by Le Doussal and Wiese in 
2009 [26] who derived the exact record statistics for a biased random walker with a 
Cauchy jump distribution (a special case of Levy flights with Levy index fi = 1). More 
recently in 2011, Wergen et al. showed that a biased random walk is useful to model 
record-breaking events in daily stock prices [11]. They were able to obtain results in some 
special limits of a biased random walker with a Gaussian jump distribution. Apart from 
these two special cases, namely the Cauchy and the Gaussian jump distribution, there 
are no other analytical results available, to our knowledge, for other jump distributions 
for a biased random walker. Recently, the record statistics for a biased random walker 
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was also studied numerically in order to quantify the contamination spread in a porous 
medium via the particle tracking simulations (27]. 

In this article we present a complete analysis of the record statistics for a 
biased random walker with arbitrary jump distributions. As we will see, the record 
statistics depends crucially on the persistence probability Q{n) [see Eq. ( TT71) below], 
the probability that the biased walker stays to the left of its initial starting position 
up to n steps. While persistence probability for various stochastic processes have been 
extensively studied in the recent past [28], it seems that for this simple biased jump 
process, it has not been systematically studied in the literature to the best of our 
knowledge. Here we provide exact results for the persistence probability Q(n) for a 
biased random walk arbitrary jump distributions [see Eq. fl67J], which subsequently 
leads to the exact record statistics for the same process. 

The rest of the paper is organized as follows. Since the paper is long with many 
detailed results, we provide in section |2]a short review on the record statistics for random 
walks both with and without bias, followed by a summary of the main results of this 
paper. Readers not interested in the details of the calculations can skip the rest of 
the paper. In section [31 we will show how to use the renewal property of the random 
walk and a generalized version of the Sparre Andersen theorem [23] to compute the 
persistence of random walks in presence of both positive and negative drift. The results 
for the persistence are interesting on their own and will be discussed in detail in section HJ 
but they will also allow us to compute the record statistics. In particular we will show 
that, in the presence of drift, the complete universality found for the record statistics 
in the unbiased case [22] breaks down and there are five different types of asymptotic 
behaviors which emerge depending on the two parameters of the model, namely the drift 
c and the index < fi < 2 characterizing the tail of the jump distribution. This record 
statistics will be discussed in detail in section [5] Later, in section |6l we will also discuss 
the extreme value statistics of the ages of the longest (section I6.ip and the shortest 
lasting records (section 16 .2 j) in each of the regimes. We will show that the asymptotic 
behavior of these quantities is also systematically different in the five regimes. Finally 
in section [7J we will conclude with some open problems. 

2. Record statistics for random walks: A short review and a summary of 
new results 

In this section, we provide a short review on the record statistics of a one dimensional 
random walk model, with and without external drift. This will also serve to set up our 
notations for the rest of the paper. At the end of this section, we summarize the main 
new results obtained in this work. 

Let us first start with the driftless case. Consider a sequence of random variables 
{xo = 0,xi,X2, ■ ■ ■ ,x n } where x m represents the position of a discrete-time unbiased 
random walker at step m. The walker starts at the origin and its position evolves via 
the Markov rule x m = x m -\ + r\ m , where r\ m represents the stochastic jump at the 
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m-ih step. The jump variables r/ m 's are i.i.d. random variables, each drawn from the 
common probability distribution function (pdf) f(rj), normalized to unity. The pdf 
f{rf) is continuous and symmetric with zero mean. Let f{k) = f{rf) e lkri drj denote 
the Fourier transform of the jump distribution. We will henceforth focus on jump 
distributions f{rj) whose Fourier transform has the following small k behavior 

f{k) = l-(l ll \k\y+... (4) 

where < \i < 2 and Z M represents a typical length scale associated with the jump. 
The exponent < fi < 2 dictates the large \r]\ tail of f{rj). For jump densities with 
a finite second moment a 2 = rf f{rf) dr], such as Gaussian, exponential, uniform 
etc, one evidently has /i = 2 and l 2 = a/y/2. In contrast, < fi < 2 corresponds 
to jump densities with fat tails firf) ~ M -1- ^ as \rj\ — > 00. A typical example is 
f{k) = exp[— \k\^] where /i = 2 corresponds to the Gaussian jump distribution, while 
< fi < 2 corresponds to Levy flights (for reviews on these jump processes see [29], [30] ) . 

A quantity that will play a crucial role later is P n {x) which denotes the probability 
density of the position of the symmetric random walk at step n. Using the Markov rule 
in Eq. (j2J), it is easy to see that P n (x) satisfies the recursion relation 

/oo 
P n ^(x')f(x-x')dx' , (5) 
•00 

starting from Po(x) = S(x). This recurrence relation can be trivially solved by taking 
Fourier transform and using the convolution structure. Inverting the Fourier transform, 
one gets 

r°° dk r a in 

Pn(x)= / ^ /(*) e~^. (6) 

In the limit of large n, the small k behavior of f(k) dominates the integral on the right 
hand side (rhs) of Eq. Substituting the small k behavior from Eq. (jlj), one easily 
finds that for < ji < 2, typically x ~ l^n 1 ^ and P n (x) approaches the scaling form [29] 

W-^(^). where ^) = £|e-l^e-^. (7) 
For < yU < 2, the scaling function C^{y) decays as a power law for large \y\ [29] 

CM — > rfer » where A » = - sm(W2) r(i + n). (8) 



A(?/) = :- - 2 . (9) 



In particular, for ji = 1, the scaling function C\{y) is precisely the Cauchy density itself 

1 1 

7r 1 + y' z 

In contrast, for /i = 2, the central limit theorem holds, x ~ a n 1 ^ 2 , and P n (x) approaches 
a Gaussian scaling form 

1 



P n (x) ^ ^j- 2 C 2 (-^) , where C 2 (y) = -j= exp(-y 2 /2) • (10) 
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From the sequence of symmetric random variables representing the position of a 
discrete-time unbiased random walker, we next construct a new sequence of random 
variables {y = 0, y u y 2 , . . . , y n } where 

y m = x m + cm implying y m = y m ^ + c + r] m , (11) 

where ?7 m 's are symmetric i.i.d. jump variables each drawn from the pdf f(rj). Clearly, 
y m then represents the position of a discrete-time random walker at step m in presence 
of a constant bias c. 

In this paper, we are interested in the record statistics of this biased sequence {yo = 
0, yx, y 2 , . . . , y n }- A record happens at step m if y m > max(y = 0, y h y 2 , ■ ■ ■ , y m -i), i.e., 
if the position of the biased walker y m at step m is bigger than all previous positions, 
with the convention that the initial position yo = is counted as a record. Let R n 
denote the number of records up to step n. Clearly, R n is a random variable and we 
denote its distribution by 

P(R, n) = Proba. [R n = R] . (12) 

We would like to compute the asymptotic properties of this record number distribution 
P(R, n) for large n, for arbitrary drift c and for arbitrary symmetric and continuous 
jump density f(rj) whose Fourier transform f{k) has the small k behavior as in Eq. (Tj0) 
with the index < \i < 2. 

In absence of a drift, i.e., for c = 0, the distribution P(R,n) was computed exactly 
in Ref. [22] for all R and n, using a renewal property of the record process. Amazingly, 
the distribution was found to be completely universal, i.e., independent of the jump 
distribution f{rj) (as long as it is symmetric and continuous) for all R and n [22]. In 
particular, for large n, it was shown that P(R, n) has a scaling form [22] 

P(i2,„)« _L g (4=) , (13) 



where the universal scaling function 



g (x) = —= exp(-x 2 /4) , for x > (14) 



is a half-Gaussian. Consequently, the mean and the variance of the number of records 
grows asymptotically as [22] 

(R n ) « 4= raV2 ' ( R D ~ ( R n) 2 « 2 ( 1 - -) n ■ (15) 

'7T \ 7T / 



The renewal property of the record process derived originally for the unbiased 
random walker in Ref. [22] was then generalized to the case with a nonzero drift c 
in Ref. [26J. In particular, the authors of Ref. [26] studied in detail the special case of 
the Cauchy jump distribution /cauch y (?7) = l/[ 7r (l + V 2 )} [which belongs to the fi — 1 
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family of jump densities in Eq. (JJ])] and found that the mean number of records (R n ) 
grows algebraically with n for large n with an exponent that depends continuously on 
c [26] 

1 11 

( R n) ~ tvi i at \\ ng(c) ' where ^( C ) = o + - arctan(c) . (16) 
1(1 + v\c) ) In 

In addition, the asymptotic distribution P(R, n) for large n was found [26] to have 
a scaling distribution, P(R,n) ~ n~ e ^ g c (i?n _e ^) with a nontrivial scaling function 
g c (x) which reduces, for c = 0, to the half-Gaussian in Eq. ( JT4l) . 

For jump densities with a finite second moment a 2 and in presence of a nonzero 
positive drift c > 0, the mean number of records (R n ) was analysed in Ref. [11] and 
was found to grow linearly with n for large n, (Rn) ~ 02(c) n where the prefactor 02(c) 
was computed approximately for the Gaussian jump distribution. However, an exact 
expression of the prefactor for arbitrary jump densities with a finite a 2 is missing. In 
addition, these results were then applied [11] to analyse the record statistics of stock 
prices from the Standard and Poors 500. The distribution of the record number P(R, n) 
for large n has not been studied for jump densities with a finite second moment. 

In this paper, we present detailed exact results for the asymptotic record number 
distribution P(R, n) for large n, for arbitrary drift c (both positive and negative) and 
for arbitrary symmetric and continuous jump densities f(r]) with Fourier transform f(k) 
having a small k behavior as in Eq. (|4]) parametrized by the exponent < fi < 2. We 
find a variety of rather rich behaviors for P(R, n) depending on the value of c and the 
exponent ji. On the strip (c, < fi < 2) (see Fig. [Q, we find five distinct regimes: (I) 
when < fi < 1 with c arbitrary (II) when fi = 1 and c arbitrary (III) when 1 < fi < 2 
and c > (IV) when /1 = 2 and c > and (V) when 1 < fi < 2 and c < 0. In these five 
regimes the record statistics behave differently, resulting in different asymptotic forms 
for the record number distribution P(R, n). The line fi = 1 (regime II above) is a critical 
line on which the record statistics exhibits marginal behavior. These five regimes are 
summarized in the phase diagram in the (c, < fi < 2) strip in Fig. [TJ 

As we will see later, a quantity that plays a crucial role in the study of record 
statistics is the persistence Q(n) which denotes the probability that the process y m in 
Eq. f lTTj) stays below its initial value y up to step n, i.e., 

Q(n) = Proba. [i/i < y , for alH = 1, 2, . . . , n] . (17) 

Due to the translational invariance of the process, Q(n) does not depend on y . The 
persistence probability has been studied quite extensively in recent years in a variety of 
theoretical and experimental systems [28] . We will see that even for the simple stochastic 
process y m representing the position of a discrete-time random walker in presence of 
a drift, the persistence Q(n) has a rather rich asymptotic behavior depending on the 
parameters [i and c. Hence, even though here our main interest is in the record statistics, 
we include the results for the persistence Q(n) as a byproduct. 

We also analyse the statistics of waiting times between individual record events. In 
particular we are interested in the expected ages of the longest and the shortest lasting 
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Figure 1. Phase diagram in the (c, < fi < 2) strip depicting 5 regimes: (I) 
< /X < 1 and c arbitrary (II) the line /i = 1 and c arbitrary (III) 1 < /j, < 2 and 
c > (IV) the semi- infinite line yu = 2 and c > and (V) 1 < ji < 2 and c < 0. The 
persistence Q(n), the record number distribution P(R,n) and the mean ages of the 
longest and the shortest lasting record exhibit different asymptotic behaviors in these 
5 regimes (see text). 



records. The age of the longest lasting record is defined as: 

/max.n = max (/i,/ 2 , •••,//?) , (18) 

where Zj is the length of the time interval between the i-th and the (i + l)-th record. 
Similarly one defines the age of the shortest lasting record as 

lmm,n = mm(h,l 2 , -;Ir) ■ (19) 

In [22], the mean values of / max ,n and Z m i n ,n were computed exactly for the symmetric 
random walk with arbitrary jump distribution. It was found that [22] for large n 

(U) ~C n, (20) 

where Co ~ 0.626508... is a universal constant independent of the jump distribution. 
Interestingly, the same constant Cq also appears in other related problems [331 El]- In 
contrast, the shortest record exhibits different behavior for large n [22] 

(^min,n 

)~V^. (21) 

In this paper we generalize these results to the case of a biased random walk and as 
in the case of record number distribution, we find five different asymptotic behaviors 
depending on c and /i. 
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Summary of the new results: Let us then summarize the main new results in 
this paper for the asymptotic behavior of the persistence Q(n), the record number 
distribution P(R, n) and the extremal ages of records in the 5 regimes in the (c, fi) strip 
mentioned above. 

Regime I (0 < fi < 1 and c arbitrary): In this regime, we find that the persistence 
Q(n) decays algebraically for large n 

Q(n) « , (22) 



where the prefactor Bj depends on the details of the jump distribution f(r]) and the 
drift c and can be computed explicitly [see Eq. ( 178 j) ]. The mean record number up to 
n steps grows asymptotically for large n as 

(R^^AjV^. (23) 

While the growth exponent 1/2 is universal, i.e. independent of c and the precise form 
of the jump distribution fijf), the prefactor A\ depends on c and on the details of 
the density f{rj). In addition, the two prefactors A\ and B\ are related simply via 
Bi = 2/(ttAi). We find the following exact expression for the prefactor A\ 



2 

Ai = —= exp 

'TT 



1 f°° dk / f(k) sm(kc) 

— arctan 



7T 



k \1 — f(k) cos (fee) 



(24) 



In the scaling limit when n — > oo and R — > oo, but with the ratio Rj ' \fn fixed, we find 
that the distribution P(R, n) approaches the scaling form 

2 / 2 R \ 1 

P(R, n)& — g I — I , where g (x) = —= exp(-x 2 /4) . (25) 
A\Ji\ n \AiWirnJ Jtx 



Averaging over R evidently reproduces the result in Eq. (|23|) . Thus, the record number, 
rescaled by the nonuniversal scale factor R — > R/A\, approaches asymptotically the same 
universal half-Gaussian scaling distribution as in the driftless case c = in Eq. f|T4|) . 

The statistics of the longest lasting record is completely unaffected by the drift c 
in this regime. For the mean value (/ m ax,n) we find that 

~ Ci n, (26) 

where the same constant C7 = C ~ 0.626508... was also found in the unbiased case [see 
Eq. (1201)]. The age of the shortest lasting record is given by 

( W> ~ A v 7 ^, (27) 

with a prefactor D\ = B\. Therefore, in contrast to (Z max ,n), (^min.n) slightly differs from 
the unbiased case and has a prefactor that depends non-trivially on c. 
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Regime II (the line \x = 1 and c arbitrary): On this line, we find that the persistence 
Q(n) decays algebraically for large n but with an exponent that depends continuously 
on c 

OH - j| > (28) 

where the exponent < 0(c) < 1 is given in Eq. (Tl6|) . In this sense the behavior is 
marginal. The pref actor Bu can be computed exactly [see Eq. (1851) ] . The mean record 
number also grows marginally for large n 

W M rIFSwr™ • (29) 

where the prefactor An = 1/ [T[l — 0(c)] Bu]. The record number distribution exhibits 
an asymptotic scaling form 

p(fl -" )ss ^ 9c (^) ' (30) 

where one can obtain a formal exact expression (11091) and explicit tails of the scaling 
function g c (x) which also exhibits marginal behavior, i.e., depends continuously on c. 

Like in regime I we find that the mean age of the longest lasting record grows 
linearly in n, but this time with a non-trivial c dependent prefactor. We find that 

(/max.n) ~ Cjl 71 , (31) 

where Cu is given in Eq. (j!47p . The mean age of the shortest lasting record is more 
strongly affected by the drift. Here we find that (/ m in,n) grows algebraically with n with 
an exponent which depends continuously on c: 

(Un,n>~ Al^ (G) , (32) 

with D\\ = Bn as in Eq. (1281) and 9 (c) as defined in Eq. ([16]) . 

Regime III (1 < fi < 2 and c > 0): In this regime, the persistence Q(n) decays for 
large n as 

Q(n) » ?f , (33) 

71^ 

where the prefactor Bm depends on the details of the jump distribution and can be 
computed [see Eq. (1901) ]. The mean number of records grows linearly with increasing n 

{R n )na^c)n, (34) 

where the prefactor a^(c) can be computed explicitly [see Eq. f )115p ]. The record number 
distribution P(R,n), for large n, behaves as 

m n) - -r^rr V » ( R ~ a ^ n ) , (35) 
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where the scaling function V^{u) can be computed exactly and it has a non-Gaussian 
form with highly asymmetric tails 

Vfj,(u) ~ | | ^ 1 as w-^-oo (36) 
« ci uP-^)/ 2 ^- 1 ) exp [-c 2 n^^" 1 )] as u -> oo , (37) 

where v4 M is the same constant as in Eq. (jS|) and the constants ci and c 2 are given 
explicitly by 

Cl = [2n^-l)^B,y/^Y 1/2 , (38) 
c 2 = (l-l/ / i)( / i J B,)^- 1 ) , (39) 

where 

2cos(/i7r/2) 

Thus, in this regime, while the mean record number grows linearly with n, the 
fluctuations around the mean are anomalous ~ n 1 ^ and described by a non-Gaussian 
distribution. 

Also the extremal ages of records have an interesting behavior in this regime. In 
particular we find that the average age of the longest lasting record grows like 

Cmax.n) ~ Cm TIM , (41) 

where the constant Cm can be computed explicitly [see Eq. (I150p ]. On the other hand 
and in contrast to the results of regime I and II, the mean age of the shortest lasting 
record converges to a finite value: 

(U,«) ~ An = 1 - a M (c) , (42) 

which thus depends continuously on a The strongly different n dependence of (/ max ,n) 
and (l m in, n ) in the regime I and in the regime III is a consequence of the fact that while 
in regime I the asymptotic behavior is dominated by the fluctuations, in regime III the 
effect of the drift is stronger in the large n limit. 

Regime IV ( the semi- infinite line \l = 2 and c > 0): On this semi-infinite line the 
variance a 2 of the jump pdf is finite. This leads to an exponential tail of the persistence 
Q(n) for large n. More precisely we show that 

Q(n)^^eM-(c 2 /2a 2 )n], (43) 

where the nonuniversal prefactor Biy can be computed exactly [see Eq. (1961)]. We also 
show that the mean and the variance of the record number both grow linearly for large n 

(R n )^a 2 (c)n and (R 2 n ) - (R n ) 2 w 6 2 (c) n , (44) 

where the amplitudes a 2 (c) and fe 2 (c) are nonuniversal and depend on the details of the 
jump distribution f{rj). We provide exact expressions for these amplitudes respectively 
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in Eqs. ( I125P and (11281) as well as in Appendix B The distribution of the record number 
P(R, n) approaches a Gaussian form asymptotically for large n 



P(R,n) pa exp 

a/27T 6 2 (c) n 



2b 2 (c)n 



(R-a 2 (c) nf 



(45) 



Thus, in this regime, the mean record number grows linearly with n with normal 
Gaussian fluctuations ~ n 1 ^ 2 around the mean. 

It is interesting to see that the asymptotic behavior of (/ max ,n) i n regime IV 
is qualitatively different from regime III. Here we find that (/ max ,n) grows only 
logarithmically with n for n — > 00: 

(U,n) ~ Civ Inn , (46) 

with an n independent constant Civ — Like in regime III, the average age of the 
shortest lasting record approaches a (different) constant value depending on c: 

(lmin,n) ~ Av = 1 - 02(c) , (47) 
which depends continuously on c. 

Regime V (1 < /1 < 2 and c < 0): In this case, the walker predominantly moves 
towards the negative axis due to the drift. Consequently, the events where the walker 
crosses the origin from the negative to the positive side become extremely rare. As a 
result, with a finite probability the walker stays forever on the negative side. Thus, the 
persistence Q(n) approaches a constant for large n 

Q(n) -> a^c) . (48) 

Similarly, the occurrence of the records (with positive record values) are also rare. 
Subsequently, we find that the mean record number also approaches a constant for 
large n 

(R n ) -> -i- , (49) 

where the constant a M (c) is given by 

a^c) = a^(\c\) for 1 < (i < 2 , (50) 
= a 2 (\c\) for (1 = 2, (51) 

where a M (c) and 02(c) are precisely the amplitude of the linear growth of the mean 
record number respectively in regime III and IV [respectively in Eqs. (134j) and ( f4*4"l) ]. 
An explicit expression for a M (c) is given in Eq. f llOip . The record number distribution 
P(R, n) also approaches a steady state, i.e., n-independent distribution as n — > 00. This 
distribution has a purely geometric form 

P(R, n^oo) = a M (c) [1 - ^(c)]*" 1 . (52) 
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Physically this result is easy to understand because for c < and [i > 1, the walker 
typically moves away from the origin on the negative side with very rare and occasional 
excursions to the positive side caused by rare large jumps. As a result, the occurrence 
of a record is like a Poisson process which eventually leads to a geometric distribution 
as in Eq. f )52l) . 

In this regime the statistics of the longest and the shortest lasting records are 
particularly simple. Since the record number is finite, the longest lasting record will 
grow linearly in n: 



with the same c dependent constant (c) as in Eq. f HHj) . Here, the main contributions 
to these quantities come from trajectories that never cross the origin and stay negative 
for all n. 

The five regimes in the (c, < fi < 2) strip are depicted in Fig. [TJ As mentioned 
above, the line \i — 1 is a special 'critical' line with marginal exponents that depend 
continuously on the drift c. It is not difficult to understand physically why /i = 1 plays 
a special role. Indeed, writing y n = x n + cn where x n represents a symmetric random 
walk, we see that the two terms x n and cn compete with each other for large n. Since 
x n ~ n 1 ^ for < jj, < 2 [see Eq. ©], it is clear that for < /i < 1, the term x n 
dominates over the drift and the presence of a nonzero drift only leads to subleading 
asymptotic effect. In contrast, for \l > 1, the drift term starts to play an important role 
in governing the asymptotic record statistics. In the region 1 < fi < 2 and c > (regime 
III), while the mean record number increases linearly with n due to the dominance 
of the drift term, the typical fluctuation around the mean is still dominated by the 
x n ~ n 1 ^ term [see Eq. (135j) ]. However when /i = 2 and c > (regime IV), the drift 
term completely dominates over the x n term leading to Gaussian fluctuations around 
the mean. This competition between x n and cn thus leads to (i) a 'phase transition' 
in the asymptotic behavior of record statistics of y n at the critical value fi = 1 and 
(ii) an anomalous region with non-Gaussian fluctuations around the mean in the region 
1 < fi < 2 and c > 0. 

3. Record Number distribution via renewal property and the generalized 
Sparre Andersen theorem 

The idea of using the renewal property of random walks to compute the distribution 
of record number was first used in Ref. [22] for symmetric random walks and was 
subsequently generalized to biased random walks [26]. We briefly summarize below 
the main idea. 

Consider the random sequence {y , y\, y 2 , ■ ■ ■ , } representing the successive positions 
of a discrete-time biased random walker evolving via Eq. ( flTT) . starting from an arbitrary 



Cmax,n) ~ C*V n 5 Cy — 1 . 

For the shortest lasting record we find a similar behavior: 



(53) 




(54) 
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initial position y . Consider first the persistence Q(n) defined in Eq. (1171) . Let us also 
define 

F(n) = Proba. \y x < y , y 2 < yo, ■ ■ ■ , y n -\ < yo, Vn > Vo] (55) 

which denotes the probability that the walker crosses its initial position y from below 
for the first time at step n. Clearly 

F{n) = Q(n - 1) - Q{n) . (56) 

It is also useful to define the generating functions 

oo oo 

Q{z) = Y,Q(n)z n , F{z) = Y,F{n)z n . (57) 

n=0 n=l 

Using the relation in Eq. (156]) it follows that 

F(z) = 1 - (1 - z)Q(z) . (58) 

Consider now any realization of the sequence {yo = 0, y±, y 2 , ■ ■ ■ , y n } up to n steps 
and let R n be the number of records in this realization. Let I = {/i, l 2 , ■ ■ ■ , Ir} denote 
the time intervals between successive records in this sequence (see Fig. [2]). Clearly 
denotes the age of the i-th record, i.e., the time up to which the i-th record survives. 
The last record, i.e. the R-th record, stays a record till step n. Let P(l, R\n) denote the 
joint distribution of the ages and the number of records up to step n. Using the two 
probabilities Q(n) and F(n) defined earlier and the fact that the successive intervals 
between records are statistically independent due to the Markov nature of the process, 
it follows immediately that 

P(l, R\n) = F(h)F(l 2 ) . . . F{l R ) Q(l R ) 5 yR N , (59) 

where the Kronecker delta enforces the global constraint that the sum of the time 
intervals equals n. The fact that the last record, i.e. the _R-th record, is still surviving 
as a record at step n indicates that the distribution Q(Ir) of Ir is different from the 
preceding ones. It is easy to check that P(l, R\n) is normalized to unity when summed 
over / and R. The record number distribution P(R,n) = J^f-P(/,-R|^) is just the 
marginal of the joint distribution when one sums over the interval lengths. Due to the 
presence of the delta function, this sum is most easily carried out by considering the 
generating function with respect to n. Multiplying Eq. ( 1591) by z n and summing over I 
and n, one arrives at the fundamental relation 

1 Q(z) , (60) 

n=0 

where we used the relation in Eq. ( |58i) . Note that, by definition, R < (n + 1), i.e. 
P(R,n) = if n < R — 1. Hence, the sum in Eq. (|60|) actually runs from n = R — 1 
to oo. 



n) z 



F{z) 



n — i _ 



Q(z)= l-(l-z)Q(z) 
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Figure 2. A typical realization of the biased random walk sequence {yo = 
0,yi,2/2) ••• ,Vn} 01 n steps with R records. Each record is represented by a filled 
circle. The set I2, ■ ■ • , Ir-i} represents the time intervals between the successive 
records and Ir is the age of the last record which is still a record till step n. 



Thus the basic object is the generating function Q(z). Once this is determined, 
one can, at least in principle, compute other quantities such as the statistics of records 
or their ages using the fundamental renewal equation f )60|) . Fortunately, there exists a 
beautiful combinatorial identity first derived by Sparre Andersen [23] that allows one to 
compute Q(z) 



Q{z) = ^Q{n) z n = exp 



n=0 



n=l 



(61) 



where p{n) = Proba. [y n < 0]. Using the relation y n = x n + cn where x n represents the 
symmetric random walk at step n in Eq. one gets, p(n) = Proba. [x n < —cn]. Then, 
using the pdf P n (x) of the symmetric walk x n at step n in Eq. ([6]), one gets 



p{n) = Proba. [x n < —cn] 



P„ (x) dx 



P„ (x) dx 



where, in obtaining the last equality we used the symmetry P n (x) 
Substituting this expression of p(n) in Eq. (IB"T]) gives 



Q(z) 



n=0 



exp 



E 

n=l 



2" 

n 



P n (x) dx 



(62) 

Pn(-X). 

(63) 



Eq. fl63|) . with P n (x) given by Eq. (jSJ), determines Q(z) in terms of the Fourier 
transform f(k) of the jump distribution fijf). Subsequently Eq. ( 1601) then determines, 
in principle, the record number distribution P(R,n). In the driftless case c = 0, great 
simplification occurs, since by symmetry f °° P n (x)dx = 1/2. This gives, from Eq. f )63l) . 
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Q(z) = 1/a/1 — z. This is completely universal as all the dependence on the jump 
distribution /(r/) drops out. Subsequently, Eq. ( 160]) provides, for c = 0, the universal 
result for the record number distribution 



Y,P{R,n)z n = ± V ) , (64) 

which, when inverted, yields (22] for large n the scaling behavior in Eq. f JT3|) with the 
scaling function given by the half-Gaussian form in Eq. ( Fl4|) . 

However, in presence of a nonzero bias c, extraction of the precise large n behavior 
of P(R,n) from the set of equations (}6"0l . (163]) and ([6]) is more complicated. For the 
special case of the Cauchy distribution, this was performed in Ref. [26] which showed 
nontrivial behavior. The rest of this paper is devoted precisely to this technical task 
of extracting the large n behavior of P(R,n) for a general jump distribution f(rf) and 
we will see that a variety of rather rich asymptotic behavior emerges depending on the 
value of the drift c and the exponent fi characterizing the small k behavior of f(k) in 
Eq. flU). 

Before finishing this section, let us remark that from Eq. ( |60|) one can also compute 
the generating functions of the moments of the number of records. For example, 
multiplying Eq. (160]) by R, summing over R and using the identity X]^Lo na;n_1 = 
1/(1 — x) 2 we get for the first moment 

Y(Rn) z n = = — . (65) 

Similarly, multiplying Eq. ( 160]) by R 2 and summing over R one gets for the second 
moment 

fVjQ^-P-'jOM. (66) 
ti (1-zfQ^z) 

We will use these two results later in Section IVB. 



4. Asymptotic behavior of persistence Q(n) for large n 

The persistence Q(n), i.e. the probability that the process y n stays below its initial value 
y up to step n and its generating function Q(z) is the key ingredient to determine the 
record number distribution P(R, n) via Eq. (|6"U]) . Apart from its key role as an input 
for the record statistics, the persistence Q(n) for this process is, by itself, an interesting 
quantity to study. We will see in this section that even for the simple stochastic process 
y n , representing the position of a discrete-time random walker in presence of a drift, the 
persistence Q(n) has a rather rich asymptotic behavior depending on the parameters ji 
and c. Before getting into the details of the derivation, it is useful to summarize these 
asymptotic results. We find that for large n, the persistence Q(n) has the following 
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asymptotic tails depending on fi > and c 

Q(n) ~ Bin~ l l 2 for < \i < 1 and c arbitrary (regime I) , 
~ B\\ n~ e ^ for /i = 1 and c arbitrary (regime II) , 
~ £>m for 1 < /j, < 2 and c > (regime III) , 
~ Bivn~ 3 / 2 exp[— (c 2 /2a 2 ) n\ for fi = 2 and c > (regime IV) , 
~ a M (c) for 1 < ji < 2 and c < (regime V) , 

(67) 

where the prefactors Bj, Bn, B m , _B IV can be explicitly computed. In regime V, a M (c) 
is a constant independent of n that can also be computed explicitly [see Eq. (110 II) and 



Appendix B for a 2 (/i)]. The exponent 6(c) depends continuously on c and is given in 
Eq. ( 1161) [see also Eq. (IHUl) ]. In Fig. [3] these results are confirmed numerically for the 
regimes I-IV. 

To derive these asymptotic behaviors of Q(n) for large n, we start with the key 
result in Eq. ( 1631) . Using Cauchy's inversion formula in the complex z plane one can 



write 



Q(n) 



J Co n=l •Jen 



(6? 



where the contour Co encircles the origin and is free of any singularity of Q(z) (see 
Fig. H]). Let z* denote the singularity of Q(z) on the real axis closest to the origin. Then, 
one can deform the contour Co to C\ (see Fig. |4j such that the vertical part of C\ is 
located just left of z* and the circular part has radius r\. By taking the r 1 — > oo limit, 
it follows from Eq. fl68l) that for large n, the contribution from the circular part vanishes 
exponentially. Thus for large n, the leading contribution comes from the vertical part 
of Ci, i.e the imaginary axis located just left of z* . Next we make a change of variable 
z = e~ sn and define 

oo 

q(s) = Q(z = e~ s ) = Q(ji) e~ sn = exp [W c »] , (69) 

n=0 

00 p— sn roo 

where W c ,„(s) = / P n {x) dx . (70) 



n 

n=l 



ci I 



Using this expression in the integrand in Eq. ( |68l) and retaining only the contribution 
from the vertical part of the contour C\ for large n, we get 



Q(n)« / —e sn exp[W c ,»(s)} , (71) 

where W c ^(s) is given in Eq. flTOl and s* = — ln(z*) is the singularity of q(s) = 
exp[W c>ll (s)} on the real axis closest to s — 0. Identifying the integral on the rhs 
of Eq. (!TT|) as a standard Bromwich integral in the complex s plane, we see that 
for large n, Q(n) is essentially given by the inverse Laplace transform of the function 




Figure 3. Numerical simulations of the persistence Q (n), i.e. the probability that a 
random walker with a bias c stays below its initial position up to step n. We considered 
4 different Levy-stable jump distributions characterized respectively by the Levy index 
/.t = 0.5, f , f .5 and /i = 2 (in the last case it is just Gaussian jump distribution). In all 
cases, we had a constant positive bias c = 1 and the data were obtained by averaging 
over 10 7 samples. For comparison, we also present the result for the unbiased case 
(c = 0) with a Gaussian jump distribution (the top curve). The thin dashed lines give 
our analytical predictions from Eq. (|67[) with fitted prefactors B\, Bu, Bm and Biy. 
For the fx = 1 case we used 9 (c = 1) « 0.7498.... 



q(s) = exp[W Cjtl (s)]. To make further progress, we need to first identify the position 
of the singularity s* of W CiJJ (s) and then analyse the dominant contribution in the 
Bromwich integral coming from the neighborhood of s* for large n. We see below that 
the singular behavior of W C)/1 (s) as a function of s depends on the parameters c and 
fi > and there are essentially 5 regimes in the (c, < /x < 2) strip as shown in Fig. [TJ 
Below we discuss these regimes separately. 

4-1. Regime I: < fi < 1 and c arbitrary 

To analyse the leading singularity of W CjIX (s) as a function of s in this regime, it is 
first convenient to use the normalization condition P n (x)dx = 1 and the symmetry 
P n (x) = P n (—x) to rewrite 

poo i pen 

/ P n {x)dx = -- / P n {x)dx. (72) 

Jen JO 
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Figure 4. The contour Co in the complex z plane can be deformed to the contour 
C\. In the large n limit, the dominant contribution to the Cauchy integral in Eq. (|68[) 
comes from the vertical part of C\ . 



Substituting this in Eq. ( J70l) gives 



i °" p — sn pen 

WcM = ~2 ln i 1 - e " s ) ~12— I. p «( x ) dx ■ ( 73 ) 

n=l 







Now, as s — > 0, the sum in Eq. ( 1?3l) converges to a constant for < /i < 1 



n /n 

71=1 >/U 



P n (x)dx. (74) 



To see this, let us see how the integral J™ P n {x) dx behaves for large n. For large 
n, we can use the scaling form for P n (x) in Eq. (j7j). One finds that J Q cn P n (x)dx — > 



(i-i/m) 



J CTl £p(y) dy as n — >• oo. For < // < 1, clearly this integral decreases leading to 
the convergence of the series in Eq. (]74| . Thus, the leading singularity of W Ctfl (s) occurs 
near s = s* = where it behaves as 

W c ^(s) « ln(s) - 5 . (75) 
Substituting this result in Eq. ( 1701) gives 

S-S>0 y/S 



n=0 



We now substitute this singular behavior of the integrand in Eq. ( 17ll after setting 
s* = and perform the standard Bromwich integral (one can use the fact that the 
inverse Laplace transform LT~^ n [s~ l ' 2 \ = l/y/wn) 

Q(n) ► A , (77) 
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By 



-S () 



7T 



exp 



-E 



n In 

n=l J ° 



P„ (x) dx 



(7? 



4-2. Regime II: \i — 1 and c arbitrary 



The case \i — 1 is rather special and marginal as we demonstrate now. Consider the sum 
W Ci i(s) in Eq. ( l70l) . In this case, it follows from Eq. (J^J) that P n {x) — > (l/n)£x(x/n) as 
n — > oo, where = 1 / [7r ( 1 + y 2 )] for all y and hence is integrable. Thus the integral 

P n (x)dx converges to a constant for large n 



/oo 
C 1 (y)dy=l-9(c) 



where 



9(c) 



Ci(y)dy — — H — arctan(c) 

2 7T 



(79) 



(80) 



Hence, the n-th term of the sum in W c ,i(s) behaves, for large n, as T n — > (1— 0(c)) e~ sn / n. 
Consequently, the sum W Ct i(s) = J2 n >i has a singularity at s = s* = 0. The leading 
asymptotic behavior near this singularity reads 

W cA (s) ► -(1 - 6(c)) \n(s) - 70 , (81) 

where 70 is a constant that depends on the details of P n (x), in particular on the difference 
between P n (x) and its large n asymptotic form (\/n)Ci(x/n) for finite n 



To 



E 

n=l 



0(c) 



P n (x) cte 



(82) 



Using this result on the right hand side (rhs) of Eq. (|70|) gives 



-70 



(83) 



Substituting this result in the Bromwich contour in Eq. ( JTTT) (after setting s* = 0) and 
performing the Bromwich integral gives 



Q(n) 



n— >oo 72 



0(c) 



where 



-To 



r[i - g(c)] and 6{c) = 2 



TV 



arctan(c) , 



34) 
35) 



and 70 in given in Eq. ( 1821 . 

Thus, for \i — 1, the persistence Q(n) decays algebraically for large n but with an 
exponent 9(c) that depends continuously on the drift c. In this sense the line /x = 1 is 
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marginally critical. The exponent 9(c) in Eq. (1851) increases continuously with c from 
9(c -> -oo) = to 9(c -»■ oo) = 1. 

Let us end this subsection with the following remark on the special case of pure 
Cauchy jump distribution, /cauchyC 7 ?) = + V 2 )}- As mentioned before, the record 

statistics for this case was studied in detail in Ref. [26]. For the Cauchy case, it is well 
known that P n (x) = (l/n)f Ca _ uchy (x/n) = (l/n)Ci(x/n) for all n. As a result, it follows 
from Eq. fl82l) that the constant 70 = in this case. However, in the general \i — 1 case 
(not necessarily the Cauchy case), 70 is generically nonzero. Thus, while the persistence 
exponent 9(c) = 1/2 + - arctan(c) is universal for all jump densities belonging to the 
/i = 1 case, the amplitude Bu is nonuniversal and depends on the details of the jump 
density. 



4-3. Regime III: 1 < yU < 2 and c > 

To analyse the singular behavior of the sum W c ^(s) in Eq. flTOj) in this regime, we 
consider the n-th term of the sum T n = (e~ sn /n) J°° P n (x)dx and substitute, for large 
n, the scaling behavior of P n (x) in Eq. (J7J). This gives T n pa (e~ sn /n) J!°(i-i/ M ) C IJ ,(y)dy. 
For 1 < < 2, the lower limit of the integral in T n becomes large as n — > 00 and we can 
use the tail behavior in Eq. (jSJ) to estimate, T n pa (A^/ fj,c' 1 )e~ sn /n M for large n. Hence 
the sum, W c ^(s) = Y^n=i T n clearly converges to a constant W Ct ^(0) as s — > 0. For 
small s, one can replace the sum by an integral and estimate exactly the first singular 
correction to this constant. This gives 



B^s 



(86) 



A»T(2 



where the constant B^ 
from Eq. (jSJ) and simplify, one finds B^ 
also that from the definition in (!70l) 



/x)/[/i(/i — l)c M ]. Using the exact expression of A^ 



-l/[2cos(/x7r/2)] > as in Eq. (SQl). Note 



q(0) = exp[W C)/ ,(0)] = exp 



00 1 



n=l 71 J 



P n (x) dx 



(87) 



Substituting the small s behavior from Eq. (!86|) in Eq. (!70|) gives 

q(s)— ->q(0)[l-B fl s^ 1 + ...} . (88) 

Substituting this singular behavior of q(s) = exp[W cll (s)] in the Bromwich integral in 
Eq. (I7ip (upon setting s* = 0) and performing the integral by standard method provides 
the following large n power law tail for Q(n) 

B\u 



Q(n) 



> 

n—>oo 



where the prefactor Bm is given by 

(»-l)B, (n-1) 



B 



111 



T(2 - n)d* 



q(0) 



2cos( / u7r/2)r(2-//)c' i 



exp 



00 1 
^— ' n 

n=l 



P„ (x) dx 



(89) 



(90) 
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4-4- Regime IV: \x = 2 and c > 

In this regime, the leading singularity s* of W C)ll (s) occurs not at s = 0, but at 
s = s* = —si where s\ = c 2 /2a 2 . To see this, let us again consider the n-th term 
of the sum W CjfM (s), i.e. T n = (e~ sn /n) P n (x)dx. For large n, P n (x) now has the 
Gaussian scaling form in Eq. (llOjl due to the central limit theorem. Substituting this 
Gaussian form and carrying out the integration one gets, 

e~ sn ( c \ 2 f°° 2 

T n — y erfc I — — yfn I , where eric(y) = —= / e~ x dx . (91) 

2n \ay/2 J y/n J y 

Using the asymptotic behavior erfc(y) ~ e~ y2 fyypK for large y, one finds that 

a e _ ( s+s i) n c 2 
T n > —= — , where Sl = — . (92) 

Consequently, the sum W c ^(s) = ^2 n>1 T n actually, while perfectly analytic near s — 0, 
has a singularity near s = s* = — s±. Close to this singular value, by taking the limit 
s + c 2 /2a 2 —> whereby replacing the sum by an integral over n, one finds the following 
leading singular behavior of W Ctfl (s) near s = — Si 

a 

«->~ai C 



W c ^{s) y Wcni-S!) -V2- , (93) 



where W c ^{— si) is just a constant. Substituting this leading singular behavior on the 
rhs of Eq. ( 1701) gives 

q(s) ► e ^M(-i) U-V2- t/s~TTi + ...] . (94) 



c 



We set s* = — Si in the Bromwich contour in Eq. (I71jl . substitute the singular behavior 
of g(s) in Eq. f )94p and perform the Bromwich integral to get 



Q{n) >^e- s ^ n where Sl = ^- (95) 

n— >oo n I 2o~ 



and the prefactor 



Biy = = — = — = exp 

cv27r cy2n 



e sin poo 



■n 



Pr, (x) dx 



(96) 



Thus, contrary to regimes I, II and III, here the persistence Q(n) has a leading 
exponential tail (modulated by a power law n~ 3 ^ 2 ). 

4-5. Regime V: 1 < \i < 2 and c < 

In this regime c = —\c\ < and /i > 1. It is convenient, using the normalization 
condition P n (x)dx = 1, to first reexpress the sum W c ^(s) in Eq. flTOj) as 



1 e -sn 



W C ,> 1 (S)=J2 / Pn(x)dx = J2 



n J-\c\n 1 n 

n=l n=l 



'OO 

/ P n (x)dx 
\c\n 



(97) 
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Performing the sum, and using the definition of W Cjfl (s) in Eq. (J70J) one gets 

W CilM (s) = - In (1 - e~ s ) - W\ 4 M . (98) 

For n > 1, Wj c | )At (0) is a constant as was demonstrated in the previous two subsections. 
Hence, one gets from Eq. ( 198]) . the leading singular behavior for small s 



W c^)^-Hs)-W Kfl (0) (99) 



which yields, via Eq. (J7D]) 



q(s) — . ^~ W ^ , (100) 

s->0 s 



Thus, in this regime, the leading singularity of q(s) occurs at s = s* = 0. Setting s* = 
and the result fllOOp in the Bromwich integral in Eq. (1711) gives 



Q(n) > a M (c) = exp[-Wjc| jA( (0)] = exp 



oo 1 
n=l 



oo 

P n (x) <ix 

|c|n 



'101] 



The fact that the persistence Q(n) approaches to a constant for large n in this regime 
can be understood physically because for c < and fi > 1, a finite fraction of trajectories 
escape to — oo as n — > oo. 

5. Asymptotic Record Number distribution P(R,n) for large n 

In this section, we analyse the asymptotic large n properties of the mean record number 
(R n ) and its full distribution P(R,n) for arbitrary c by analysing the set of equations 
©, (1601) . (1631) and (1651) with arbitrary jump distribution f{rf). Consider first the mean 
record number. As in Section IV, we invert Eq. ( 1651) by using the Cauchy inversion 
formula, deform the contour (as in Fig. HJ), keep only the vertical part of the contour C\ 
for large n and finally make the substitution z = e~ s to obtain the following Bromwich 
formula 

ns ' +io ° ds .„ 1 



{R ^I (!-,-.).,-(.) ' (1 ° 2) 



S*— IOO 



where q(s) is given in Eqs. ( |69|) and ( 1701) and its small s properties have already been 
analysed in section IV in different regimes in the (c, < fi < 2) strip. As in section IV, 
s* denotes the singularity of q(s) on the real line in the complex plane that is closest to 
the origin at s — 0. 

Similarly, the record number distribution is obtained by inverting Eq. ( l60l) in the 

same way 

P(R jn )^ / ^e sn q(s) [l-Cl-e-)^)]*" 1 . (103) 

In this section, we use the already derived results for q(s) in Section IV and analyse 
the asymptotic behavior of (R n ) and P(R, n) respectively in Eqs. (1 1Q2[) and ( I103[) in 
different regimes of the (c, < fi < 2) strip and on the critical line \i — 1. 
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Figure 5. a): Rescaled mean record number (R n )/^/n for a Levy-stable distribution 
with Levy index /i = 1/2 and different series length n = 10 2 ,10 3 ,10 4 and 10 5 . For 
each n the average was performed over 10 3 samples. For n 3> 1 the results collapse and 
agree with the predicted analytical behavior for A\ (c) in Eq. (|104|) . b): Rescaled 
distribution Aiy/wh P(R,n) /2 as a function of 2R/ Ai^/wn of the record number R n 
after n steps for a random walk with a Levy-stable jump distribution of Levy index 
(j, = 1/2, n = 10 3 and different values of the drift c = —1,0,0.1,1 and 10. We also 
plotted the asymptotic analytical result go(x) given in Eq. (|2"5j) . All curves collapse 
nicely. In regime I, the record number has a half-gaussian distribution. 



5.1. Regime I: < // < 1 and c arbitrary 

Let us first consider the asymptotic behavior of the mean number of records (R n ) for 



large n in this regime. Consider the Bromwich integral in Eq. 



For large n, this 



integral can be shown to be dominated by the small s region of the integrand. Taking 
s — > limit in the integrand, substituting the result f lTo]) on the rhs of Eq. (11021) . and 
performing the Bromwich integral we get the leading asymptotic behavior for large n 



(R n ) w Aiy/n , where Ai = —= e s ° = —= 



7T 



7T 



exp 



n=l 



n 



P„ (x) dx 



(104) 



Comparing this to the amplitude of persistence in Eq. ( !7H|) we see that the two prefactors 
are related simply via B\ = 2/(itAi). The prefactor A\ can further be expressed explicitly 
in terms of the Fourier transform of the jump distribution f(k) as in Eq. (124j) . This 
is shown in Appendix A where we also compute the asymptotic behavior of A\ for 



In Fig. [5] a) we compare this result for {Rn) to numerical 

fn agree nicely with our analytical 



large |c| [see Eq. ( 1A.9I) 

simulations. The numerical results for fi>l, {Rn) 
values for Aj (c). 

Next we turn to P(R, n) in the limit of large n. To extract the scaling behavior of 
P(R, n) from Eq. (I103p . we substitute on the rhs the small s behavior of q(s) from Eq. 
( 1761) and use the notation e~ s ° = (2/y/Tr)Ai. The appropriate scaling limit is clearly 
R — )■ oo, s — > but keeping the product s/sR fixed. Taking this limit in Eq. ( 1103P 
gives, 



P(R, n) 



ds 
2jti 



AiJk~s 



exp 



A l 



7T 



yfsR 



(105) 
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c 



Figure 6. ln(i?„) / In n as a function of the drift c for the Cauchy distribution with Levy 
index /i = 1 and for different values of n = 10 2 , 10 3 , 10 4 and 10 5 . For each n and c, the 
average was performed over 10 3 samples. The results from the numerical simulations 
collapse and agree with the predicted analytical behavior of ln(i? n )/lnn = 9 (c) and 
9 (c) = i + ^arctan (c) as in Eq. ([85]>. 

One can simply evaluate the Bromwich integral by using the identity, LT~J^ n [e~ bRy ^ j 'y/s\ = 
e -b 2 R 2 /&n i This leads to the asymptotic result announced in Eq. (1251) in the scaling 
limit n — > oo, R — > oo with the ratio R/y/n fixed. In Fig. [5]b) we computed numerically 
the rescaled distribution Aw/wn P(R,n) /2 as a function of IRj ' Awprm and compared 
it with go(x) Eq. (125]) . The figure confirms that in regime I, the record number has a 
half- Gaussian distribution with a width that depends non-trivially on the drift c and 
the Levy-index fi. 

In summary, for < \i < 1, the drift is not strong enough to change the \fn growth 
of the mean record number. The presence of drift just modifies the prefactor of the \fn 
growth. Similarly, the distribution of the record number in Eq. (125]) in presence of a 
drift, when appropriately scaled, remains unchanged from the universal half-Gaussian 
form in the driftless case. 

5.2. Regime II: \x = 1 and c arbitrary 

As mentioned in the introduction, on the critical line \i — 1, the record statistics 
was investigated in detail in Ref. [26] for the special case of Cauchy jump distribution 
/cauchyl??) = VKC^ + ^ or a g enera l jump distribution with /i — 1 (not necessarily 
of the Cauchy form), the record statistics has a very similar mathematical structure 
that can be derived from the general framework developed in this paper. 

Let us first consider the growth of the mean record number (R n ) in Eq. (110211 . 
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Substituting the small s behavior of q(s) from Eq. ( I83p and performing the Bromwich 
integral upon setting s* = we get for large n 

(Rn) « rn Al \ ( v . n e ^ where A u = . (106) 
T(l + 9(c)) 

Note that 70 is a distribution dependent constant while the exponent 0(c) = 
P_ Ci(y)dy — 1/2 + - arctan(c) is universal. In Fig. (JBJ) this exponent is plotted and 
compared with numerical simulations of random walks with a Cauchy jump distribution 

(/*=!)■ 

Turning now to the distribution P(R,n) in Eq. f 1 1 U 3 j) . as before, we substitute the 
small s expansion of q(s) from Eq. (1831) . It turns out that the appropriate scaling limit 
for P(R, n) is n — > 00, R — >■ 00 but keeping the ratio R/n 9 ^ fixed. To see this, we first 
set s* = 0, set R large but fixed, and keep the leading terms for small s to get 

/+«°o J„ 1 
77- e Sn -TZ^ exp l-e~^ s 6 ^ R] . (107) 
•too ^ m S 

Rescaling sn 4 s and keeping the scaled variable R/n 6 ^ fixed gives the asymptotic 
scaling distribution 

P (R, n) « -r^r 9c (-r^) , (108) 



A lin 6 ^ * c \A n n d ( c ) 

where the scaling function g c {u), which depends continuously on c, is given by the formal 
Bromwich integral 

to 8 9(c)-l e .-uSl°) with M > Q; (1Q9) 
-too 2m 

where we recall that < 9(c) < 1. 

One can easily extract the tail behavior of the scaling function g c (u) by analysing 
the integral in Eq. f 1 1 9 j) . For instance, when u — > 0, g c (u) approaches a constant 

/+ioo j -1 1 

, ioo ^^'^-^^i^^nT^M' (110) 

The integral in Eq. f II 1 j) can be performed by wrapping the contour around the branch 
cut on the negative real s axis. 

In the opposite limit, when u — > 00, the integral in Eq. (I109p can be performed 
using the standard steepest descent method. Skipping details and using the shorthand 
notation 9 = 9(c) we get 

g c (u -+ 00) « [2tt(1 - 0) fl^/M] " 1/2 u -(i-2«)/2(i-«) exp [-(1-0) ^/d^) u i/(i-*)] . 

(Ill) 

Thus the distribution has a non-Gaussian tail. The function g c (u) can be expressed in 
terms of the one-sided Levy distribution, which was discussed for instance in Ref . [32] . In 
some particular cases, the Bromwich integral in Eq. (11091) can be evaluated explicitly. 
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1 2 3 4 5 



R/A|, n e(c) 

Figure 7. Rescaled plot of An n e ^ P(R, n) as a function of R/Aun 6 ^ for /i = 1 and 
c = -1/V3 = -0.57735, and 9{c) = 1/3 (re gime II). These data have been obtained 
for n = 10 5 and averaged over 10 5 samples. The dotted line corresponds to our exact 
result in Eq. (fTT2|) . 



For rational values of 0(c), g c {u) can be expressed as a finite sum of hypergeometric 
functions. A very special case corresponds to c = —l/y/3 where one has 9 = 1/3, such 
that 

9c=-i,A u ) = 32/3 Ai (3573) > u > • ( 112 ) 
where Ai(x) is the Airy function. Its asymptotic behaviors are then given by 

9<=-iW5(«) ~ Vr(2/3) , u -» (113) 



^W»)~S"- V4 «p(-^). d«) 

which agree with the general analysis presented above flll0[ imp . In Fig. [7] we show 
a plot of the rescaled probability Aun 8 ^ P(R,n) as a function of R/Ann e ^ computed 
numerically for c = — l/v3, which agrees reasonably well with our exact analytical 
result in Eq. ffTT2l . 




Figure 8. Rescaled distribution a fl (c)n 1 / tl P{R^n) of the record number R n after 
n steps for a random walk with a Levy-stable jump distribution of Levy index 
fj, = 1.5. The data are plotted as a function of the shifted and scaled variable 
u = (R — a^c) n)/(a Al (c) n 1 ^). The different curves correspond to different values 
of n = 10 3 , 10 4 , 10 5 and 10 6 and for a drift c = 1. They were obtained by averaging 
over 10 6 samples. For n = 10 5 and n = 10 6 the numerical results were binned for 
technical reasons. We also plotted our analytical results for the scaling function V M (it) 
given by Eq. (|121j) . While for smaller values of n, there is still a significant difference 
between the simulations and our analytical result, it converges to the behavior in 
Eq. (|! 21j) when n increases. 



5.3. Regime III: 1 < /i < 2 and c> 



We first compute the asymptotic growth of the mean number of records in this regime. 
Substituting the leading singular behavior of q(s) from Eq. flHHj) on the rhs of Eq. d!02[) 
and performing the Bromwich integral gives 



(R n ) ~ a M (c) n where 



1 



9(0) 



exp 



^ n 

n=l 



P n (x) dx 



(115) 



Note that we used above the expression of q(0) in Eq. flH7|) . We have checked numerically 
this linear growth and in Fig. [10] the bottom curve shows a plot of (R n ) /n as a function 
of c, although we have not tried to evaluate a M (c) numerically. 

We next consider the distribution P(R,n) in Eq. (I103p . We substitute the small 
s behavior of q(s) from Eq. (1HH1) on the rhs of Eq. (11031) . set s* = 0, R large and keep 
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only leading small s terms to get 



r +tco ds 

P(R,n)&q(0) / exp [s (q(0)R - n) + BM)Rs*\ . (116) 

J-ioo 2m 

Next we set 

R = afl (c)n + a^c) n 1/fl u , (117) 

where a M (c) = 1/^(0) and take the limit R — > oo, n — > oo but keeping the scaled variable 
u above fixed. We substitute Eq. (I117P on the rhs of Eq. (j!16p . Keeping only the two 
leading terms for large n and fixed u gives 



/ +tco ds 
exp l-sn^^u + B,n s^} . 
-foe 2m 



'MS) 



Note that for fixed u, both terms inside the exponential are of the same order. In fact, 
the scaling in Eq. (j!17p is chosen so as to make the two leading terms precisely of the 
same order for large n. Rescaling s by n 1//M , i.e., sn 1 ^ — > s and using a^(c) = l/q(0) 
reduces Eq. (11181) to a nicer scaling form announced in Eq. (|35p 



P(R,n) 



V.Au), where u - 



and the scaling function V^(tt) is formally given by the Bromwich integral 



119) 



ds 
2tH 



-us+Bu 



(120) 



where the constant > is given in Eq. ( 1401) . 

Interestingly, the same scaling function V^(w) also appeared in Ref. [3T] in the 
context of the partition function of the zero range process on a ring. The asymptotic tails 
of the function V^iu) were analysed in great detail in j3l] (see Eqs. (78)- (83) and Fig. 5 
in Ref. [31] and note that in [31], the index \i was denoted by 7 — 1). We do not repeat 
the computations here, but just quote the results. It was found that V^u) has highly 
asymmetric tails. For u — > —00, it decays as a power law, V^(u) — > K [J \u\~ IJ '~ 1 where the 
prefactor = B^T^l+fi) sm[ii(fi+l)]/n. Using our expression B^ = — 1/(2 cos(/i7r/2)) 
from Eq. ( )40|) . it is easy to show that = where the constant is defined in Eq. 
(JHJ). This leads to Eq. ( 136]) . In contrast, when u — > 00, V^(u) has a faster than Gaussian 
tail as described precisely in Eq. ( 1371) . To plot this scaling function, a convenient real 
space representation can be used from Ref. [31]. Replacing 7 — 1 by fi in Eq. (84) of 
Ref. [31] and using B^ = — 1/2 cos(/i7r/2), we obtain 



71 



dye 



-2/72 



cos 



tan(/i7r/2) y^ + yu 



'121) 



We compared this result for a Levy index of /x = 1.5 to numerical simulations in Fig. |HJ 
Even though the convergence of the numerically obtained distributions is slow, it is clear 
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Figure 9. Rescaled distribution a fl (c)n 1 / tl P{R^n) of the record number R n after 
n = 10 4 steps for a random walk with a Levy-stable jump distribution with different 
Levy indices /j, = 1.25,1.5,1.75 and \i = 2. The data are plotted as a function of 
the shifted and scaled variable u = (R — a Al (c) n)/(a A1 (c) n 1 ^). For all these data, the 
value of the drift is c = 1 and they have been obtained by averaging over 10 6 samples. 
The figure shows that for fi — > 2 this rescaled distribution approaches the Gaussian 
form given in Eq. (|134l) . 



that the asymptotic distribution (u) is approached for n — > oo. In Fig. [9] we plotted 
numerical simulations of the rescaled record number distribution for different values of 
/i. One finds both numerically and by taking the limit in Eq. (I134p that, for \i — \ 2, 
this rescaled distribution approaches a Gaussian form (see regime IV). 

To summarize, in this regime the mean record number increases linearly with 
increasing n, but the typical fluctuations around the mean are anomalously large 
of Ofo 1 ^) (superdiffusive) as described in Eq. ( 11171) . In addition, the probability 
distribution of these typical fluctuations around the mean are described by a highly 
non- Gaussian form described precisely in Eq. (11191) . 



5.4- Regime IV: = 2 and c > 

In this regime, as explained in section IV.C, q(s) = exp[W c ^(s)} in Eqs. fl69|) and flTD]) 
is analytic at s = 0. This can be seen by expanding the sum W c ^{s) in Eq. (1701) in a 
Taylor series in s 

W Cjlx {s) = J2 d m s m , where d m = ^-L- ^ n m ~ l / P n {x) dx . (122) 

m=0 m ' n=l -* cn 
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Figure 10. Numerical simulations of (R n )/n for random walks with a Gaussian 
(with variance a = 1), an exponential [with parameter b = 1, see its definition below 
Eq. (|126j) ]. both regime IV, and a Levy-stable jump distribution with /j = 1.5, in regime 
III, with positive drift c > 0. For each distribution we show data for n = 10 4 which 
were obtained by averaging over 10 4 samples. For the Gaussian and the exponential 
distribution we also plotted a numerical evaluation of our exact formula for a,2 (c) using 
Eq. (|126|) for the Gaussian case and Eq. (|B.12|) for the exponential case. Those curves 
agree perfectly with the numerical simulations. 



The coefficient d mj for each m, is finite as the sum over n is convergent since the integral 
fen Pri(. x ) decreases with n faster than exponentially for large n (see section IV. C), 
as long as /x = 2 and c > 0. Consequently, for small s, q(s) also has a Taylor series 
expansion 

q(s) = q(0) + q'(0) s + V(0)s 2 + . . . (123) 



q(0) + q'(0)s+-q"(0)s 2 + ... 

Let us start with the asymptotic behavior of the mean record number (R n ) in 
Eq. (11021) . Once again, the dominant contribution to the integral in Eq. ( 1102}) for large 
n comes from the small s region. Taking the s — > limit in the integrand and using 
the small s expansion in Eq. (11231) . keeping only the leading terms and performing the 
Bromwich integral term by term one gets for large n 



(Rn) 



where 



ds 
2tu 



g(0)s 2 



(1 



g(0) 



0(s s 



a 2 [c)n + k 2 [c) 



a 2 (c) 



and «2(c) = [1 



g(0) 

7(0)/?(0)]/g(0). 



exp 



oo 1 

^ n 

71=1 



P n (x) dx 



0(l/n) 
(124) 

(125) 
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{2Tia 2 Y 1 / 2 e-^/ 2a \ we have 
1/2 e -x 2 /2a 2 n an j [j ence 02(c) in Eq. (11251) is given by the explicit 

(126) 



For example, for a Gaussian jump distribution f(j]) 
P n {x) = (2vma 2 
formula 

00 

c\/n 



a 2 (c) = exp 



a 



V2 



— — erfc 
^ 2n 

n=l 

For instance, for c = 1 and a = 1, one gets a 2 (c = 1) = 0.800543 .... Another example 
is the exponential jump distribution f(rj) = (2 6) _1 exp(— \x\/b). In this case, one can 
also compute (see the Appendix B) the constant a 2 (c) = A where A is given by the 
solution of the transcendental equation exp(— Ac/6) = 1 — A 2 . For example, for c = 1, 
6=1, one gets A = 0.714556 .... For these two examples, we have confirmed the leading 
asymptotic result for the mean record number in Eq. ( 11240 with the exactly computed 
prefactors a 2 (c) (as discussed above) in our numerical simulations (see Fig. [TOl . 

In a similar way, one can also analyse Eq. (166]) for the large n behavior of the 



second moment (-R 2 ). Skipping details, we get the following leading large n behavior 



(R 2 n ) » a 2 2 (c)n 2 + p 2 {c)n + 0(l) , where p 2 (c) 



3 - q(0) - 4^ 



<f(o) L 

Consequently, the variance of the record number grows for large n as 



9(0) J 



(R 2 n ) - (R n ) 2 « b 2 {c)n where b 2 {c) 



q 2 (0) 



9(0) 



, g(o) 

' 9(0) 



(127) 



'1281 



Thus, in this regime, the mean record number grows linearly with n for large n while 
the size of typical fluctuations around this mean grows as ~ \J~ri. 

How are these typical fluctuations around the mean distributed? To answer this, 
we need to analyse P(R, n) in Eq. ( 11031) in the scaling limit where both n and R are 
large, but the ratio (R — a 2 (c)n) / \fn is fixed. To proceed, we set s* = and substitute 
the small s expansion of q(s) in Eq. (11231) on the rhs of Eq. (11031) . take R large but 
fixed to get 



ds 
2m 



exp [s (q(0)R 



n) 



(1/2) b 2 (c)q 3 (0)Rs 2 ] 



where 6 2 (c) is given in Eq. (11280 . Next we set 

R = a 2 (c) n + \Jb 2 (c) \fnu 



(129) 



(130) 



where a 2 (c) = 1/^(0) is given in Eq. (1 1 2 5 [) and take the scaling limit where R —> oo, 
n — )• oo but keeping the scaled variable u above fixed. Substituting R from Eq. ( 11301) 
into Eq. (11291) and keeping only the two leading terms for large n gives 



P(R,n) w q(Q) 



+ IO0 



ds 
2~Ki 



exp 



VMS) 9(0) V^s u + (1/2) 62(c) q 2 (0)n s' 



(131) 



Note that for fixed u, both terms inside the exponential are of the same order. Indeed, 
as in the section VB, the scaling in Eq. ( I130p is chosen so as to make the two leading 
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Figure 11. Plot of the cumulative distribution of record numbers P < (R,n) = 
Proba. [R n < R] as a function of the shifted and scaled variable u = 
(R — 02(c) n) /(-\/^2(c) n) for a random walk with Gaussian jump distribution (with 
a = 1) of n = 10 4 steps. The different curves correspond to different values of positive 
drift c = 1/16, 1/4, 1 and 2. For each c the data were obtained by averaging over 10 6 
samples. We compared the numerical results to the cumulative distribution of V2 (n), 
which we obtained analytically (Eq. (|134|l ). All curves collapse nicely, confirming that 
the asymptotic record number of a biased Gaussian random walk with a positive drift 
has the Gaussian distribution given by Eq. (|132l) . 



terms precisely of the same order for large n. Rescaling \/b^{c)q(0)y^n s — > s simplifies 
to 

R — a 2 (c) n 



P(R,n) 



V 2 {u) where u 



\fb 2 (c)n y / b 2 {c)n 
and the scaling function V 2 (u) is given by the Bromwich integral 

' 3 ds 



VoU] 



-us+s 2 /2 



2ni 



which can be exactly computed (since it is a Gaussian integral) to give 

1 



VM 



^7T 



exp[— u /2] . 



(132) 



(133) 



(134) 



This then proves that P(R,n) is asymptotically Gaussian as announced in Eq. (I45j) . 
Fig. m~1 confirms this result numerically. We plotted the cumulative distribution of record 
numbers P< (R, n) = Proba. [R n < R] as a function of the shifted and scaled variable 
u = (R — 02(c) n) /(\/b 2 (c) n) after n = 10 4 steps for different values of positive drift c 
and compared them to a Gaussian cdf (cumulative distribution function). All numerical 
results collapsed perfectly on the analytical curve. 
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0.5 1 1.5 2 2.5 3 3.5 4 

u = R * a 2 (|c|) 



Figure 12. Rescaled distribution a,2(\c\)P(R 1 n) of the record number R„ after 
n = 10 steps for a random walk with a Gaussian jump distribution, of variance a = 1, 
with different negative values of the drift c = — 0.01,c = —0.05,-0.1 and —0.25. 
The data are plotted as a function of the rescaled variable u = Ra 2 {\c\). For each 
value of c the data were obtained by averaging over 10 samples. We compared the 
numerical results with a simple geometric distribution. The good agreement confirms 
our analytical findings given by Eq. (|136[) . 



5. 5. Regime V: 1 < /i < 2 and c < 

In this regime, we set s* = in Eq. ( 1 1 3 f> and substitute on its rhs the small s expansion 
of q(s) from Eq. (llOOp . Keeping only leading order behavior for small s gives, for large n, 

P(R,n)^a,( C )[l-a,(c)] R ^ / — e sn - , (135) 

where the constant a M (c) = exp[— W| c | iM (0)] = exp 
Eq. (flUT]) . 

Using the fact that LT~^ n [l/ s\ = 1 gives the large n (but R fixed) behavior of 
P(R,n) 

P(R, n) > a M (c) [1 - a^c)]*- 1 . (136) 

n—>oc 

Thus, the distribution becomes independent of n for large n and has a simple geometric 
form with mean (R n ) —> l/a^(c). Comparing the expression of a^(c) as given in 
Eq. f llOip and those of a^(c) in Eq. (II 15ft and a 2 (c) in Eq. (11251) for c > 0, one 
immediately finds that a M (c) = a M (|c|) for 1 < ji < 2 while a 2 (c) = a 2 (\c\), the results 
mentioned respectively in Eqs. ( |50l) and ( 15"TT) . 

In Fig. [12] we compared Eq. (1 1 3 6 [) to numerical simulations of negatively biased 



as given m 
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Gaussian random walks with different values of c. For large n the rescaled distribution 
of u = Ra2(\c\) approaches the geometric (exponential) distribution e~ u . 



6. Extreme statistics of the age of a record 

From the previous study of the mean number of records (R n ), one deduces that the 
typical age (see Fig. [2])) of a record is given by / typ ~ n/(R n ). However, following 
Ref. [22] for the unbiased case, it turns out that the extreme ages of records do not 
share the typical behavior. In this section, we probe such atypical extremal statistics 
by considering the longest and shortest lasting records characterized by their respective 
ages (durations) I 

max,n and /min,n- 'We focus on their mean values (/max,n); (^mm,n) and 
find rather different asymptotic behaviors in the five regimes in the (c, < \i < 2) strip 
mentioned before (Fig. [T|). 



6.1. Age of the longest lasting record / max ,n 

We first consider the longest lasting record whose age / max ,n is given by (see Fig. [2]) 

/ max = max(/i,/ 2 , • • • ,Ir) ■ (137) 

The cumulative distribution J-* n (m) = Proba. (Z max ,n < m) was studied in Ref. [22], 
where an explicit formula for its generating function (GF) was obtained: 

f>M*" = jlSu » (138) 
n=o Z^=i I ) 

where F(l) = Q(l — 1) — Q(l), from which one deduces the generating function of the 
mean (Z max ,„) = X)~ =1 [l - F n (m)} 

E y n/j \_S^\ 1 YmLiQ(1) z1 

Z \'max,n/ / 



n=0 

OO v-^OO 



feii-^ i-Er^^'j (139) 

i V YZ m F{i)z l + {i-z)Y.T= m Q(iV (uo) 

where we have used that F(z) = 1 — (1 — z)Q(z) (!58l) . 

In the absence of drift, c = 0, it was shown in Ref. [22] that (l max ,n) behaves, 
for large n, linearly with n with a non trivial coefficient, independently of the jump 
distribution f{rj) 

</,„.„„> ~ C n , C = [ dy —^-±—^- = 0.626508... (141) 

Interestingly, this constant C appears also in the study of the longest excursion of 
Brownian motion [331 EI]. Note that to obtain the large n behavior of (/ m ax.n) from 
Eq. (11391) one has to analyse the above formula fl 1 3 9 [) in the limit z — > 1. We will see 
that in this limit the above sum over m is dominated by the large values of m, which 
thus depends crucially on the large m behavior of the persistence probability Q(m). 
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Consequently (/ max ,n) behaves quite differently in the five regimes in the (c, < /x < 2) 
strip in Fig. [1] and are summarized as follows: 



(^max,n) ~ ti for < fi < 1 and c arbitrary (regime I) , 
~ n for fi = 1 and c .arbitrary (regime II) , 
~ n~ L for 1 < ji < 2 and c > (regime III) , 
~ In n for /i = 2 and c > (regime IV) , 
~ n for 1 < ji < 2 and c < (regime V) . 



(142) 



In the following we will discuss the behavior of (/ m ax,n) separately for the five regimes. 



6.1.1. Regime I: < fi < 1, c arbitrary: In this regime, we remind that Q(m) behaves, 
for large m, as 

Bi s B i 



Q(m) 



F(m) 



(143) 



'm 2m 3 / 2 ' 

where Bi is given in Eq. (1781) . Setting z = e~ s we are interested in the limit s — > in 
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Figure 13. Plot of (Z maXi „)/n in the different regimes I, II and V: the points are the 
results of our numerical simulations. For regime II (/i = 1), wc present two curves, one 
with a negative drift (c = —1) (the second curve from top) and one with a positive drift 
(c = 1) (the bottom curve). These data indicate that in all these cases (i max .„) oc n, for 
large n, with an amplitude which agree quite well with our analytical results, which 
are represented in solid line for each of these cases and corresponds to the formula 
given in Eq. $M fT47lfT58)) . 
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the formula in Eq. (11391) where one can replace F(m) and Q(m) by their asymptotic 
behaviors 

oo oo 1 v^oo /-3/2 p -sl , „ V^oo 7-1/2 „-sl 

n=0 m=l v 1 i—di—ra 

where we have used q(s) ~ \pnB\j \fs when s — > (!76| [7H|) . In the limit s — )■ 0, the 
discrete sums over I and m can be replaced by integrals and one finds that the right 
hand side in Eq. (I144p behaves like 1/s 2 when s — > with a prefactor which we can 
compute to obtain the large n behavior of (Z max ,n) as 

(/max,n) ~ Ci 71 , Cj = / dy , f oo _ 3/2 _ x = C , (145) 

where C is given above (I14ip and where the last equality is simply obtained by 
performing an integration by part in the integral over x in the denominator. In Fig. [131 
we have plotted the results of our numerical estimate of (£ max ,n) (obtained by averaging 
over 10 4 different realizations of random walks) for \x = 0.5 and two different values of 
c = ±1.0. This plot shows that (l ma . x ,n)/n saturates rather quickly to the constant Co, 
independently of c, in agreement with Eq. (j!45p . 

Thus in this regime the large n behavior of (Z maXi n) is unaffected by the presence of 
the drift c. This result could have been anticipated as /! max ,n can be considered as the 
longest excursion between two consecutive zeros of a renewal process with a persistence 
exponent 1/2. This quantity was studied in Ref. [31] and its average was computed, 
yielding the large n behavior obtained in Eq. (j!45p . 



6.1.2. Regime II: fj, = 1 and c arbitrary: In this regime, we recall that the persistence 
probability Q(m) behaves algebraically for large m with an exponent 9(c) which depends 
continuously on c 

B 11 

( 3( m ) Wc) ' ^( c ) = o + ~ arctan ( c ) > ( 146 ) 

where the amplitude Bu is given in Eq. (|85|) . Here again we can use the result obtained 
in Ref. [M] for the longest excursion between consecutive zeros of a renewal process 
with a persistence exponent 9(c) to obtain 

f°° 1 
(U ^ Cn "' Cu = J dy i + y^ e y^dxx-^e-^ (147) 
which depends continuously on c and is independent of the non-universal amplitude 
Bu (1146ft . In Fig. IT4l we show a comparison of Cfi obtained numerically (the squares 
symbols) and from our exact formula (solid line), which shows a very good agreement 
between both. 



6.1.3. Regime III: 1 < ji < 2 and c > 0: In this regime the persistence probability 
Q(m) behaves for large m as 

Q(m) ~ ^ , (148) 
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Figure 14. Plot of C\\ as a function of c. The red squares correspond to numerical 
data while the solid line corresponds to our analytical result in Eq. (|147|) together 
with Eq. dUBJ). 



where the amplitude Bm is given in Eq. 
Ref. EU one obtains that 



(lmax,n) ~ C m ™ 



i/m 



( 190]) . Using again the results obtained in 



(149) 



where, however, the amplitude Cm was not given in Ref. [34]. A careful analysis of the 
above formula ( I139D allows to obtain the amplitude Cm as 



C 



in 



T(l - I//*) 



(EL 



— sm , 

vr V 2 



r(A*) 



(150) 



which diverges as Cm ~ (7r(/x — l))" 1 when /i — y 1 and vanishes as Cm ~ a/7t(2 — /x)/2 
when /i — j- 2. In Fig. [T5]we show a plot of our numerical data for (/ max ,n) (averaged again 
over 10 4 different realizations) for different values of \i = 1.4, 1.5, 1.7, 1.9 and for a fixed 
value of the drift c = 5.0. The solid lines indicate the corresponding exact asymptotic 
behaviors in Eq. (I149[ 1150]) : the agreement between the two is quite good although the 
convergence to the asymptotic behavior gets slower as \i decreases to 1. 



6.1.4- Regime IV: \i — 2 and c > 0: In this case the persistence Q(m) behaves quite 
differently as it vanishes exponentially for large m as 



Q(n) 



B 



n 



IV 

3/2 



-sin 



where s\ 



2a 2 



(151) 
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Figure 15. Plot, in a log-log scale, of (l maxn ) as a function of n in regime III: the 
different curves correspond to different values of fj, = 1.4, 1.5, 1.7, 1.9 with a fixed value 
of c = 5.0. The solid line are the exact results given in Eqs (|149[ I150j) . without any 
fitting parameter. 



where the amplitude Biy is given in Eq. f )96|) . This case was not analyzed in Ref. [3 
From Eq. (JT39J) one has in this case 

Therefore in the limit when s — > one can estimate the leading behavior of the sum 
over m as 

m 
s 



V(/ max ,n)e- sn , (153) 

n=0 

where m* is such that 



Q{m*) ~ sg(0) . (154) 

2 

From the asymptotic behavior above (I15ip one finds that m* ~ — ^ his so that finally 

2a 2 

(/max,n) ~ ClV hi 71 , Cry = —5- , (155) 

which is in sharp contrast with the algebraic growth obtained above in Eq. (I149p for 
1 < fi < 2 and c > 0. In Fig. [JB] we show a plot of (/ max ,n) as a function of lnra: 
the straight line suggests indeed a logarithmic growth, in agreement with our analytic 
result (I155p . However, a more precise comparison with this exact asymptotic result, as 
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Figure 16. Plot of (l maxn ) as a function of Inn in the regime IV: here \x = 2 
and the two curves correspond to c = 1 and c = 1.5 (c = 1 in both cases). The 
two curves suggest a logarithmic growth, as expected from Eq. ()155|) . Inset: Plot 
of (?max,n) — 2 In n where 2 Inn is the exact asymptotic result from Eq. (|155|) and 
2cr 2 /c 2 = 2. This plot suggests rather strong corrections oc In Inn to the leading 
logarithmic growth of (Z ma x,n)- 



shown in the inset of Fig. [161 suggests that the leading corrections are proportional to 
In In n, and hence quite strong. 



6.1.5. Regime V: 1 < \i < 2 and c < 0; In this case the persistence probability Q(m) 
tends asymptotically to a constant (1101 j) : 



Q{m) > a M (c) = exp[-W\ c \ jfM (0)] = exp 



°° -I POO 

V- / P n {x)dx 



(156) 



In addition from (1561) one has that Q(m) — a^(c) oc n^ 1 so that F(m) oc m~ fl for large 
m. Therefore, the terms entering into the sum in Eq. (11 39[) are given, to leading order 
when 1 — z — e~ s — > and large m (which are terms which give the leading contribution 
to this sum over m) 



Therefore this yields 

(/max,n) ~ Cy 71 , Cy = 1 . 



q(0) 



■e ~ sm = e~ sm 



(157) 



(158) 
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This result, which is corroborated by our numerical simulations (see Fig. IT3"j) . can be 
physically understood as in this regime where c < and \i > 1 the number of records is 
finite and these records typically occur during the first steps of the random walks, where 
the walker might stay positive for a short while before it escapes to negative values when 
n — > oo, and no record happens any more. 

6.2. Age of shortest lasting record / m i njn 

We now consider the shortest lasting record whose age l m m,n is given by (see Fig. [2]) 

^min,n = min(/i,/ 2 , • ■ ■ • (159) 

Note that, given that the final incomplete interval Ir is taken into consideration above, 
/ min n can be zero: this happens when a record has been broken at the last step, such 
that In = 0. 

The cumulative distribution Q n (m) = Proba. (/ m i njn > m) was studied in Ref. [22] 
and an explicit formula was obtained for its generating function: 

p^^ i E M%* > (160) 

from which one gets the generating function of the average value (Z m m,n) as 



In the absence of drift, c = 0, it was shown in Ref. [22] that 



(161) 



(W) ~ D ^ > D = ~i=- ( 162 ) 

\/7T 



As for (/ max ,n) we will see that the behavior of (l m i n ,n), in the presence of non zero drift 
c 7^ 0, is quite different in the five different regimes discussed above. Again we start by 
giving a brief summary of our results for (Z m i n ,n): 

(^min,n) ~ V™ f° r < fi < 1 and c arbitrary (regime I) , 
~ n 1 " 9 ^ for fi = 1 and c arbitrary (regime II) , 
~ const, for 1 < /i < 2 and c > (regime III) , 
~ const, for fi = 2 and c > (regime IV) , 
~ n for 1 < /i < 2 and c < (regime V) , 

(163) 

again with 6(c) as defined in Eq. ( 1161) . In the following we discuss the behavior of 
(^mm,n) in more detail for each of the five regimes. 

6.2.1. Regime I: < fi < 1 and c arbitrary In this case the persistence probability 
decays algebraically as given in Eq. ( II 43ft and the analysis of (l m m, n ) can be obtained by 
noticing that, in the limit z — > 1, the denominator in Eq. (|16ip can be simply replaced 



Record statistics and persistence for a random walk with a drift 



12 




n 



Figure 17. Plot, on log-log scale, of (l m i n ,n) as a function of n, for different values of 
(i < 1 and c (regime I) . The points are the results of numerical simulations while solid 
lines correspond to our exact analytic result given in Eq. (|164|) . These data indicate 
that in this regime (i m m,n) y^n, although the corrections to the exact asymptotic 
behavior are clearly visible, in particular for fj, = 0.8, c = 1.0. 



by 1 while the remaining sums over / (in the numerator) and over m can be replaced by 
integrals. This yields straightforwardly 

</min,n> ~ AV^ , (164) 



L>i = B t = —= exp 

'71 




f(k) sin (kc) 



(165) 

1 — f(k) cos (kc) 

where the expression of B\ is given in Eq. ( I7H1) . In Fig. [T71 we show the results of our 
numerical simulations which are in a rather good agreement with Eq. (I164p . although 
the corrections to this exact asymptotic behavior are clearly visible, in particular for 
fi = 0.8, c = 1.0. In Fig. \TE\ we show a plot of the numerical computation of (/ m m )ri ) for 
fi = 1 and different values of c = —1, 0.5 and c = 1: these data are in good agreement 
with the power law growth in Eq. ( I164p . although we have not attempted to estimate 
numerically the pref actor Dj. 

6. 2. 2. Regime II: \x = 1 and c arbitrary In this regime where the persistence probability 
Q(m) decays algebraically as in Eq. ( 1841) . (/ m m,n) can be analyzed as in the regime I 
where in the limit z — > 1, the denominator in Eq. f II 6 1 j) can be simply replaced by 1 
while the remaining sums over I (in the numerator) and over m can be replaced by 
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Figure 18. Plot, on a log-log scale, of (l m i n ,n) as a function of n for /i = 1 and different 
values of c = —1,0.5 and c = 1. The solid line corresponds to the algebraic growth 
n i-9(c) ) from Eq (fi66]> . 



,2-0(c) 



OO />00 



integrals. This yields straightforwardly 

OO R 

V e - ffl (I ■ ) ~ — 11 

/ o \ l mm,n/ 2- 
n=l 

which yields 



cfy / dxx- 9 <- e) e- v 



B 



ii 



,2-«(c) 



r[2-%)] 



(/min,n)~Al^ (c) , D n = B n , 

where -Bn is given in Eq. (1831) and 0(c) = 1/2 + - arctan(c). 



(166) 



(167) 



6.2.3. Regime III: 1 < \i < 2 and c > In this case we write the above formula (11611) 



as 



1 / 1 \ 

Si . <Uw> (! - M )+Ei 

n=0 v ^ " m=2 



oo v-^oo 



168) 



where we have simply isolated the term m = 1 and used 1 — F(0) = (1 — z)Q(0) ( 158]) . 
Now the above sum fl!68j) . which starts with m = 2, is dominated by the large values of 
m. Because of the algebraic decay of Q(m) ~ m _At in this case ( I148P and fi > 1 in this 
regime one gets that this second term behaves like (1 — z)^~ 2 , which is then subleading, 
compared to the first term which behaves like (1 — z) -1 . Therefore one gets in this case 



(lmin,n) ~ D 



III , J^lll 



g(0) 



exp 



- / P n (x) dx 
z — ' n 

n=l 



(169) 
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Figure 19. Plot of (Z m in,rc) as a function of n for /j, = 1.5 and \i = 2 and different values 
of c > 0, therefore corresponding to regime III and IV. The solid line corresponds to 
the exact result, from Eq. (fl69l [T7D|> . 



where we have used the expression for 1/^(0) given in Eq. (11151) . In Fig. [19] we show a 
plot of the numerical computation (/ m m !n ) for [i = 1.5 and different values of c = 0.5 and 
c = 1, which is in very good agreement with Eq. ( 1 1 6 9 f) . Note that we have extracted 
the value of 1/^(0) which enters into the expression of Dm from the linear growth of 
the mean record number (R n ), according to (11151) . 



6.2.4- Regime IV: fi = 2 and c > A similar analysis can be carried out in this case, 
starting from the same formula (11681) . In this case, in the above sum (11681) . which starts 
with m = 2, one can safely put z = 1, because of the behavior of the exponential decay 
of Q(m) in this case (I15ip . Therefore one gets immediately 



(4nin,n) ~ Dfy , Djy 



1 



q(0) 



1 — cxp 



E- / P n (x)dx 
n 



n=l 



(170) 



where we have used the expression for l/q{0) given in Eq. (I125p . In Fig. [19] we show a 
plot of the numerical computation (/ m in in ) for fi = 2 and c = 1, which is good agreement 
with Eq. (I170p . Note that we have extracted the value of 1/^(0) which enters into the 
expression of Diy from the linear growth of the mean record number (R n ), according to 
Eq. (H25D. 
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Figure 20. Plot of (Z m m,ra) /ft as a function of n for different values of 1 < fj, < 2 and 
different values of c < 0, corresponding to regime V. 



6.2.5. Regime V: 1 < fi < 2 and c < In this regime where the persistence goes to 
a constant <5(m) — > a M (c), for m 3> 1, one can simply replace Q(l) by this constant 
value in the sum of the numerator in Eq. (jlfiip while the denominator can be simply 
approximated by 1 in the limit 1 — z = e~ s — > 0. This yields straightforwardly 

(£min,n) ~ <ty(c) " ■ (171) 

In Fig. 120] we show a plot of (/ m m,n)/^ which we have computed numerically for different 
values of \x = 1.7, 1.5 and fi = 2 and also for different values of the drift. These results 
are in very good agreement with our exact asymptotic result in Eq. ( 1171 jl . where the 
value of a M (c) have been extracted from the mean record number (R n ) ~ l/a M (c) f T49|) . 
This result ( jl Tilt can be easily understood by realizing that l m in,n — n if the whole 
trajectory is on the negative side, which happens with probability q; m (c) while l m m,n is 
of order 0(1) if the walker makes an excursion on the positive side. One also notices 
that, in this case, / typ ~ (/ m in,n)- 

7. Conclusion 

In this paper we considered a very simple model of a one dimensional discrete-time 
random walk in presence of a constant drift c. At each time step the particle jumps by 
a random distance c + t] where the noise r) is drawn from a continuous and symmetric 
jump distribution f(t]), characterized by a Levy index < fi < 2. The jump has a 
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Figure 21. The figure shows numerical results for the mean record number (R n ) for 
biased random walks from all five regimes. For regimes I to IV we used a positive 
bias of c = 1, in regime V we simulated a Gaussian random walk (with a = 1) with a 
negative bias of c = —0.01. For each jump distribution we averaged over 10 4 samples. 
In all these cases, as shown in detail in the previous figures, the asymptotic behavior 
agree very well with our analytical predictions (which are not shown on this figure for 
clarity). 



finite second moment for /x = 2, while for < /i < 2 the second moment diverges. For 
this discrete-time series consisting of the successive positions of the biased walker, we 
presented complete analytical studies of the persistence and the record statistics. For 
the later, we studied the mean and the full distribution of the number of records up to 
step n and also the statistics of the duration of records, in particular those for the longest 
and shortest lasting records. As a function of the two parameters c and < \x < 2, 
we found that it is necessary to distinguish between five different universal regimes, as 
shown in the basic phase diagram in Fig. [TJ In these 5 regimes, the persistence and 
the record statistics exhibit very different asymptotic behaviors that are summarized in 
Section 2 and we do not repeat them here. For instance, the growth of the mean record 
number with n in all five regimes is summarized in the simulation results in Fig. ETJ in 
complete agreement with our analytical predictions. The main conclusion is that even 
though this is a rather simple model, it exhibits very rich and varied universal behaviors 
for record statistics and persistence depending on the two parameters c and < fi < 2. 
Our results provide a simple yet nontrivial, but fully solvable model for the record 
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statistics, a subject which has gained considerable interest over the last few years. 
Our results provide one generalization of the previous results for record statistics for 
symmetric random walks (22]. However, it is important to note that this extension does 
not yet cover all possible kinds of discrete-time random walks. In principle one could 
consider more complicated asymmetries of the jump distribution. It might be interesting 
to consider a jump distribution that has different tail-exponents in the left and in the 
right tail. Also a generalization of these results to an asymmetric lattice random walk 
is still missing. In [22] a symmetric lattice random was also considered. It should be 
possible to compute the record statistics of a lattice random walk that has a higher 
probability to jump in one direction than in the other. 

It might be interesting to see if our results can be applied to financial data, similar to 
the analysis in (TTJ [12] . Daily stock data however proved not to be useful for comparison 
because the asymptotic limit is hardly achieved in the available observational data. An 
application to stock data with a higher temporal resolution however should be possible 
and might provide new insights. Such an analysis is definitely an interesting subject for 
future research. Also the distribution of records in stock prices has not been analysed 
in detail before and it would be interesting to see if such an analysis for available data 
can be fitted to our theoretical distributions. 
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Appendix A. The constant A\ 

The constant A\ in Eq. H104I) can be directly expressed in terms of f(k) as announced 
in Eq. (124"]) . To derive this, we use the explicit expression of P n {x) from Eq. ([6]) in the 
expression for A\ and integrate over x to get 




exp 




(A.l) 



Next we use the symmetry f(k) = f(—k) which leads to 




exp 



7T 



1 



f 00 dk ^ sin(kcn) ~ 

J iL^T L/WJ 

J0 n=l 



n 



(A.2) 



The sum on the rhs can be explicitly evaluated using the identity 




n=l 



(A.3) 
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which then leads to the exact expression in Eq. 

We then analyze the behavior of A\ when |c| is large and in the case where 
f(k) = exp(— l&l^ 1 ), with fi < 1. In that case one has P n (x) = n'^-^C^x/n 1 ^) for 
all n and it is easier to start from the formula given in the text in Eq. (11041) 



9 °° i pen 

A = -^e 5 ° , S = S (c) = J2 ~ / ^(x/n^dx/n 1 ^ . 



(A.4) 



Note that, given that P n (x) = P n (—x) one has Sq(c) = So(— c) and we thus present the 
analysis for c > 0. Performing the change of variable y = x/n 1 ^ in the integral above 
(IA.4j) we write 



ccn v- 



_ i pen f 

So(c) = J2~ £Mdy, (A.5) 
and take the derivative with respect to c 

OO , s 

*) = £«-^ 4, . (A.6) 

n=i Vn f J 

In this expression, one notices that c/n~ = (n/c 1 ^)" so that when c — > oo the 
discrete sum over n in Eq. ( 1A.6I) can be replaced by an integral (we recall that fi < 1 
here), which leads to 

W ~\J~c* (y^) y~ lhl d y • ( A - 7 ) 

Finally, performing the change of variable z = y~ in Eq. ( 1A.7I) yields 

cl-piJo c 2(1-/4) 



so that one gets 



2 m 

^ = — e So oc c 5 ^ , c ->• oo . (A.9) 

/7T 



This power law behavior (IA.9j) can be understood from the following scaling argument. 
We are indeed interested in the records statistics of the variables y n , with y n = 
x n + cn ( TIT]) where x n behaves for large n as x n = Ofa 1 ^). Therefore for small n, n < n* 
when c is large, y n is dominated by the drift term and n* is such that cn* ~ n* 1 ^, which 
yields 

n* ~ cw^ . (A.10) 

On the other hand, for small n, n < n* , y n is dominated by the (positive) drift and 
hence is almost deterministic which yields (R n ) ~ n, for n < n* while (-R„) ~ Ai^fn for 
n> n* . By matching these two behaviors for n = n* one obtains 

Ay ~ Vn* oc c^C^I , (A. 11) 

which yields the result obtained above ( 1A.7j) . 

Note finally that, by using S (c) = —S (—c) one obtains 

A l rsj (-c) 5(h5 , c -> -oo . (A. 12) 
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Appendix B. Computation of 02(c) = ei2(|c|), c < for exponential jump 
distribution with c < 

The expression for the amplitude 02(c) in regime V (with c < 0) and for a general 
jump distribution is given in Eq. ( 11011) . By comparing with Eq. ( 1125ft we see that 
oj2(c < 0) = 0,2 ( I c| ) where 02 ( | c| ) is the prefactor of the leading linear growth of mean 
record number in regime IV with drift positive \c\. For a general jump distribution f(rj), 
we then have 



Fourier transform of the jump distribution. Thus, in general, computing the prefactor 
02 (c) = 02 ( |c| ) explicitly is difficult for arbitrary f{rf). It can be done explicitly for 
Gaussian distribution where P n (x) = (27ina 2 )~ 1 ^ 2 exp[— x 2 /2na 2 ] itself is Gaussian and 
02(c) = a 2 ( I c| ) is then given by the formula in Eq. (I126p . In this appendix, we show 
that 02(c) = a 2 ( I c| ) can also be computed explicitly for the symmetric exponential 
distribution f{rj) = (2 6) _1 exp(— \x\/b). 

For this exponential jump distribution, the Fourier transform has the Lorentzian 
form, f{k) = l/[ii(b 2 k 2 + 1)]. One can then substitute this in the expression for P n (x) 
and eventually in Eq. (IB.lj) . After a quite convoluted computation involving contour 
integration in the complex plane, one can find 02(c) explicitly. However, as we show 
below, for the exponential case, there is an alternative simpler way to compute 02(c) 
directly (without going through the formula in Eq. ( IB. II) . 

The first observation is that 02(c) is just the limiting value of the persistence 
probability Q{n) (the probability that the walker stays below up to n steps starting 
at 0) when n — > 00 in presence of a negative drift c < 0. By symmetry, Q(n) is then 
also the probability that the walker, starting at the origin, stays above the origin up 
to n steps, but in presence of a positive drift |c| > 0. So, the idea is to compute this 
probability Q(n) directly for the exponential jump distribution and then take the limit 
n 00 to compute 02(c) = Q(n — )■ 00). 

To compute Q(n), we first define 

In (v) = Proba. that the random walker, starting at y > stays positive up to step n . 



If we can compute ^(y), then Q{n) is simply obtained by putting the starting position 
to be 0, i.e., Q{n) = g^(0). To compute q^iy), we can write a backward recurrence 
relation for q^{y) by considering the jump that happens at the first step from y to 




(B.l) 



where we recall that P n (x) = |£ \f(k)\ e 



—ikx 



and f(k) = f '(rj) e tkri dr\ is the 



(B.2) 



y'>o 




(B.3) 



Qo(y) = 1 fory > 0. 



(B.4) 
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In the limit of large n, we expect that q^iu) approaches to an n independent stationary 
value, q^iy) q + (y), that just denotes the eventual probability with which the walker 
escapes to infinity (starting from y) in presence of a positive drift |c|. Taking n — > oo 
limit on both sides of Eq. f ]B.3j) gives the integral equation for y > 



q + (y)= / q + (y')f(y+\c\-y')dy'. (B.5) 



o 



Note that this equation is valid for arbitrary jump distribution /(//). This half-space 
Wiener-Hopf type integral equation with asymmetric kernel can not be solved in general. 
However, for the special case of the exponential distribution, f(r)) = 1/ (2b) exp( — J77I/6), 
this integral equation (IB.5j) can be transformed into a differential equation using 



f"(v) = ~S(r,) + i/fo) . (B.6) 
By differentiating twice Eq. ( 1B.5|) with respect to y one then obtains [using Eq. (1B.6j) 

^0^ = -^ + (y + \c\) + ^q + (y). (B.7) 

Note that the solution q + (y) must approach to 1 as y — > 00: q + (y — > 00) = 1. This 
follows from the fact that if the particle starts at the positive infinity, it escapes to 
positive infinity with probability 1 in presence of any positive drift |c| > 0. 

Note that the differential equation (1B.7|) , though linear, is actually nonlocal in y due 



to the first term on the rhs and hence is still not completely trivial to solve. Fortunately, 
it turns out that it admits a solution of the form 

q + (y) = l-Fexp(-Xy/b), (B.8) 

where F and A are two dimensionless constants (independent of y) that are yet to 
be determined. Note that this ansatz manifestly satisfies the boundary condition 
q + {y ~~ 00 ) = 1- Substituting this ansatz in Eq. (IB. 7ft we see that indeed Eq. (IB.8P is 
a solution provided A satisfies the equation 

exp (-A \c\/b) = 1 - A 2 ; with A > . (B.9) 

The transcendental equation has a unique positive solution which then determines A 
uniquely. For example, for b/c = 1, we get using Mathematica the root A = 0.714556 . . .. 
But we still need to determine the prefactor F in the ansatz in Eq. (IB. 81) . The amplitude 
F in Eq. (IB.8j) is obtained by injecting this solution back into the integral equation (IB. 71) 



and performing the integral. Indeed, one finds that Eq. ( 1B.8j) is a solution of the integral 
equation provided 

F=l-A. (B.10) 
This then uniquely determines the solution of the integral equation (1B.7j) 



q+{y) = l-(l-X)exp(-\y/b) (B.ll) 



where A is the unique positive solution of the transcendental equation (IB.9H . 
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Noting finally that a 2 (c) = Q(n — > oo) = g + (0) gives 

a 2 (c) = a 2 (\c\)=q + (0) = X, (B.12) 

where A > is the solution of Eq. ( IB. 91) . We have checked that we indeed get exactly 
the same expression by evaluating the original general expression in Eq. (IB. II) for the 
exponential jump distribution, though this was not completely trivial to check (we do 
not give details of this check here). 
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